MECHANISM  OF  INTERMITTENT  ATOMIZATION 
FINAL  REPORT 


DTIC 

electe  ' 

JUL  0  9 1993 


A 


S.  P.  Lin 
June  1993 

U.S.  ARMY  RESEARCH  OFFICE 
DAAL03-89-0179 
CLARKSON  UNIVERSITY 


APPROVED  FOR  PUBLIC  RELEASE 


DISTRIBUTION  UNLIMITED 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  r*oort»ng  buroe*  tor  thu  collection  of  information  n  eitimatw  to  *vr»9f  '  hour  par  r«por*e.  including  the  time  for  reviewing  instruction*.  searching  ousting  oata  sources 
gathering  and  maintaining  the  data  needed,  and  co meeting  and  renewing  the  collection  of  information  Send  comments  reoardmg  thi»  burden  estimate  or  anv  other  aspen  of  t*«s 
collection  of  information,  including  suggestions  tor  reducing  this  burden  to  Washington  MeaOQuanea  Services,  Directorate  for  information  Operations  and  Heooni.  1215  Jefferson 
Oevts  Highway.  Suite  1204.  Arlington,  vA  22202-002.  andtotne  Office  o»  Management  and  Budget.  Paperwork  deduction  Project  (0704-0 i»).  Washington,  DC  20503 

1.  AGENCY  USE  ONLY  (Leave  blink)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  OATES  COVERED 

June  15,  1993  Final  Report,  10/1/89  -  6/30/93 

4.  TITLE  AND  SUBTITLE 

Mechanism  of  Intermittent  Atomization 

S.  FUNDING  NUMBERS 

DAAL03-89-K-0179 

S.  AUTHOR(S) 

S.P.  Lin 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS(ES) 

Clarkson  University 

Potsdam,  NY  13699 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

370350F 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AODRESS(ES) 

U.  S.  Army  Research  Office 

P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 

10.  SPONSORING  /  MONITORING 

AGENCY  REPORT  NUMBER 

2<*i  H-t  0'E<r 

11.  SUPPLEMENTARY  NOTES 

The  view,  opinions  and/or  findings  contained  in  this  report  are  those  of  the 
author(s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 

12».  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

A  theory  is  developed  for  the  purpose  of  enhancing  the  efficiency  of  internal 
combustion  engines,  and  to  reduce  the  air  pollution  caused  by  inefficient  combustion.  The 
theoretical  results  show  that  the  droplet  formation  from  the  leading  edge  of  the  liquid  fuel  jet 
is  caused  by  the  Taylor-Lamb  instability  during  the  decerelation  phase  of  the  intermittent 
spray.  The  mechanism  of  droplet  formation  from  the  rest  of  the  liquid-gas  interface  is  caused 
by  the  interfacia]  stress  fluctuation.  The  decerelation  of  the  liquid  jet  tends  to  reduce  the 
droplet  diameter  at  the  same  pressure.  However,  finer  sprays  can  be  produced  with  a  more 
sudden  reduction  in  the  nozzle  inlet  pressure  after  its  peak  during  the  injection  period.  An 
increase  in  the  compressibility  of  the  liquid  jet  is  shown  the  promote  the  generation  of  smaller 
droplets,  although  at  a  slower  rate.  The  compressibility  of  the  ambient  gas,  as  well  as  the 
surface  tension,  tends  to  promote  absolute  instability.  The  precise  physical  meaning  remains 
unclear.  The  absolute  instability  can  be  reduced  to  convective  instability  by  reducing  the 
surface  tension  and/or  increasing  the  compressibility  of  the  jet  fluid.  The  convective 
_ instability  is  responsible  for  the  formation  of  sprays. _ 

14.  SUBJECT  TERMS 

Intermittent  fuel  sprays,  atomization 

IS.  NUMBER  OF  PAGES 

88 

IB.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION  IB.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION 

OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED  UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

UL 

MSN  7540-01 -280-5500  Standard  form  298  (Rev  2*89 


*r«cne*«J  by  4*5'  Std  23»-iS 


m-102 


GENERAL  INSTRUCTIONS  FOR  COMPLETING  SF  298 


The  Report  Documentation  Page  (RDP)  is  used  in  announcing  and  cataloging  reports.  It  is  important 
that  this  information  be  consistent  with  the  rest  of  the  report,  particularly  the  cover  and  title  page. 
Instructions  for  filling  in  each  block  of  the  form  follow.  It  is  important  to  stay  within  the  lines  to  meet 
optical  scanning  requirements. 


Block  1.  Agency  .  Use  Only  (Leave  blank 


Block  2.  Report  Date.  Full  publication  date 
including  day,  month,  and  year,  if  available  (e.g.  1 
Jan  88).  Must  cite  at  least  the  year. 

Block  3.  Type  of  Report  and  Dates  Covered. 


State  whether  report  is  interim,  final,  etc.  If 
applicable,  enter  inclusive  report  dates  (e.g.  10 
Jun  87-30Jun88). 

Block  4.  Title  and  Subtitle.  A  title  is  taken  from 


the  part  of  the  report  that  provides  the  most 
meaningful  and  complete  information.  When  a 
report  is  prepared  in  more  than  one  volume, 
repeat  the  primary  title,  add  volume  number,  and 
include  subtitle  for  the  specific  volume.  On 
classified  documents  enter  the  title  classification 
in  parentheses. 

Block  5.  Funding  Numbers.  To  include  contract 


and  grant  numbers;  may  include  program 
element  number(s),  project  number(s),  task 
number(s),  and  work  unit  number(s).  Use  the 
following  labels: 


Contract 

Grant 

Program 

Element 


Project 

Task 

Work  Unit 
Accession  No. 


Block  6.  Author(s).  Name(s)  of  person(s) 
responsible  for  writing  the  report,  performing 
the  research,  or  credited  with  the  content  of  the 
report.  If  editor  or  compiler,  this  should  follow 
the  name(s). 

Block  7.  Performing  Organization  Name(s)  and 


Address(es).  Self-explanatory. 

Block  8.  Performing  Organization  Report 


Number.  Enter  the  unique  alphanumeric  report 
number(s)  assigned  by  the  organization 
performing  the  report. 

Block  9.  Sponsoring/Monitorinq  Agency  Name(s) 
and  Address(es).  Self-explanatory. 

Block  10.  Sponsoring/Monitorinq  Agency 
Report  Number.  (If  known) 


Block  11.  Supplementary  Notes.  Enter 


information  not  included  elsewhere  such  as: 
Prepared  in  cooperation  with  ..;  Trans,  of...;  To  be 
published  in....  When  a  report  is  revised,  include 
a  statement  whether  the  new  report  supersedes 
or  supplements  the  older  report. 


Block  12a.  Distribution/Availabilitv  Statement. 


Denotes  public  availability  or  limitations.  Cite  any 
availability  to  the  public.  Enter  additional 
limitations  or  special  markings  in  all  capitals  (e.g. 
NOFORN,  REL,  ITAR). 


See  DoDD  5230.24,  'Distribution 
Statements  on  Technical 
Documents." 

See  authorities. 

See  Handbook  NHB  2200.2. 

Leave  blank. 


DOE 

NASA 

NTIS 


Block  12b.  Distribution  Code. 


NASA- 
NTIS  • 


Leave  blank. 

Enter  DOE  distribution  categories 
from  the  Standard  Distribution  for 
Unclassified  Scientific  and  Technical 
Reports. 

Leave  blank. 

Leave  blank. 


Block  13.  Abstract.  Include  a  brief  (Maximum 
200  words)  factual  summary  of  the  most 
significant  information  contained  in  the  report. 

Block  14.  Subject  Terms.  Keywords  or  phrases 
identifying  major  subjects  in  the  report. 

Block  15.  Number  of  Pages.  Enter  the  total 


number  of  pages. 

Block  16.  Price  Code.  Enter  appropriate  price 
code  (NTIS  only). 

Blocks  17.  -19.  Security  Classifications.  Self- 


explanatory.  Enter  U.S.  Security  Classification  in 
accordance  with  U.S.  Security  Regulations  (i.e., 
UNCLASSIFIED).  If  form  contains  classified 
information,  stamp  classification  on  the  top  and 
bottom  of  the  page. 

Block  20.  Limitation  of  Abstract.  This  block  must 


be  completed  to  assign  a  limitation  to  the 
abstract.  Enter  either  UL  (unlimited)  or  SAR  (same 
as  report).  An  entry  in  this  block  is  necessary  if 
the  abstract  is  to  be  limited.  If  blank,  the  abstract 
is  assumed  to  be  unlimited. 


Standard  Form  298  Back  (Rav  2-89) 


1 


MECHANISM  OF  INTERMITTENT  ATOMIZATION 
FINAL  REPORT 

* 

S.  P.  Lin 
June  1993 

U.S.  ARMY  RESEARCH  OFFICE 
DAAL03-89-0179 
CLARKSON  UNIVERSITY 

APPROVED  FOR  PUBLIC  RELEASE 
DISTRIBUTION  UNLIMITED 


TVT7C’  QUALITY  INSPECTED  8 


FOREWORD 
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1.  Problem  Statement 


The  fundamental  mechanism  of  the  formation  of  intermittent  sprays,  which  are  encountered 
for  example  in  diesel  engines,  is  analyzed  theoretically.  The  effects  of  the  variation  of  the 
pressure  wave  form,  which  causes  the  intermittency,  on  the  characteristics  of  sprays  is 
investigated.  The  compressibility  of  the  ambient  gas  as  well  as  that  of  the  jet  liquid  are  taken 
into  account  The  speed  of  the  jet  ranges  from  subsonic  to  supersonic. 

2.  Summary  of  Research  Results 

It  is  shown  that  the  droplet  formation  from  the  leading  edge  of  the  liquid  fuel  jet  is  caused 
by  the  Taylor-Lamb  instability  associated  with  the  apparent  body  force  during  the  deceleration 
phase  of  the  intermittent  spray.  The  mechanism  of  droplet  formation  from  the  rest  of  the  liquid- 
gas  interface  remains  the  same  as  that  for  the  case  of  stationary  atomization,  i.e.  the  interfacial 
stress  fluctuation.  It  is  found  that  finer  sprays  can  be  produced  with  a  more  sudden  reduction 
in  the  nozzle  inlet  pressure  after  its  peak  during  the  injection  period.  An  increase  in  the 
compressibility  of  the  liquid  jet  is  shown  to  promote  the  generation  of  smaller  droplets,  although 
at  a  slower  rate.  The  compressibility  of  the  ambient  fluid,  as  well  as  the  surface  tension  are 
shown  to  promote  absolute  instability.  The  absolute  instability  can  be  reduced  to  the  convective 
instability  by  increasing  the  compressibility  of  the  jet  liquid.  The  convective  instability  is  the 
manifestation  of  spatially  growing  disturbances.  The  physical  significance  of  the  theoretically 
predicted  absolute  instability  remains  poorly  understood.  An  effort  by  the  present  P.I.  to 
understand  this  significance  is  currently  supported  by  a  NASA  grant.  The  viscosity  of  the 
ambient  gas  is  demonstrated  to  increase  the  droplet  diameter  in  a  spray.  These  results  are 
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published  in  journal  papers  given  in  Appendices  A-D.  Appendix  E  "Experimental  Study  of  an 
Intermittent  Spray  Using  Several  Nozzle  Pressure  Waveforms"  is  excerpted  from  the  M.S.  thesis 
of  Vincent  F.  Cook.  He  is  expected  to  write  up  his  thesis  for  journal  publication  in  the  near 
future.  The  theoretical  counter  part  of  his  thesis,  the  excerpt  of  which  is  given  in  Appendix  F, 
is  being  prepared  for  publication. 
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6.  Appendices 


A 


Reprinted  from  June  1992,  Vol.  114,  Part  2.  Journal  of  Applied  Mechanics 


(APPENPIX  A  ) 


E.  A.  Ibrahim 

Mechanical  Engineering  Department, 
Tuskegee  University, 
Tuskegee.  AL 
Mem  ASME 


S.  P.  Lin 

Department  of  Mechanical  and 
Aeronautical  Engineering. 
Clarkson  University. 
Potsdam,  NY 
Mem  ASME 


Weakly  Nonlinear  Instability  of  a 
Liquid  Jet  in  a  Viscous  Gas 

The  weakly  nonlinear  instability  of  a  viscous  liquid  jet  emanated  into  a  viscous  gas 
contained  in  a  coaxial  vertical  circular  pipe  is  investigated  as  an  in  it  tab  value  problem. 
The  linear  stability  theory  predicted  that  the  jet  may  become  unstable  either  due  to 
capillary  pinching  or  due  to  interfacial  stress  fluctuation.  The  results  cf  nonlinear 
stability  analysis  shows  no  tendency  of  supercritical  stauJity  for  both  of  the  linearly 
unstable  modes .  In  fact,  the  nonlinear  growth  rate  of  the  disturbance  is  faster  than 
the  exponential  growth  rate  of  the  linear  normal  mode  disturbance  for  the  same 
flow  parameters.  Moreover,  the  most  amplified  linear  normal  mode  disturbance 
evolves  nonlinearly  into  a  nonsinusoidal  wave  of  shorter  waielength.  t\o  nonlinear 
instability  is  found  frr  the  linearly  stable  disturbances.  Thus,  while  the  linear  theory 
is  adequate  for  the  p  edict  ion  of  the  onset  of jet  breakup ,  nonlinear  theory  is  required 
to  describe  the  outcome  of  the  jet  breakup. 


Introduction 

A  computational  investigation  of  the  stability  of  cylinderical 
jets  is  undertaken  as  an  attempt  to  understand  the  atomization 
phenomena  of  liquid  jets.  These  phenomena  occur  in  many 
applications  such  as  the  liquid  fuel  injection  in  internal  com- 
bustion  engines,  spray  coatings  of  surfaces,  forming  submicron 
solid  panicles,  and  in  various  printing  devices. 

For  many  purposes,  the  linear  theory  adequately  describes 
the  salient  features  of  the  jet  instability;  for  example,  it  ac¬ 
curately  predicts  the  small  amplitude  gr^w  ’  rate  of  a  sinu¬ 
soidal  disturbance.  However,  it  is  a  commor  observation  that 
the  wave  form  of  an  initially  sinusoidal  perturbation  becomes 
nonsinusoidal  close  to  the  point  of  drop  formation,  i.e.,  the 
point  at  which  the  surface  wave  amplitude  equals  the  radius 
of  the  jet  (in  capillary  instability).  In  order  to  explain  this 
phenomena,  much  work  has  been  done  on  the  nonlinear  cap¬ 
illary  instability  of  liquid  jets.  Wang  (1968)  investigated  the 
effect  of  a  finite  amplitude  on  a  circular  column  of  inviseid 
liquid  with  a  third-order  perturbation  approach  but  points  out 
that  this  solution  is  not  uniformly  valid  for  all  wave  numbers. 
Yuen  ( 1 968)  also  used  the  perturbation  approach  to  third  order, 
but  his  choice  for  the  third-order  term  gives  rise  to  a  secular 
term  whose  elimination  leads  to  a  small  amplitude-dependent 
shift  in  the  cutoff  wave  number.  Yuen  (1968)  concluded  that 
the  asymmetrical  development  of  an  initially  sinusoidal  wave 
is  a  nonlinear  effect  with  generation  of  higher  harmonics  as 
well  as  feedback  into  the  fundamental.  Nayfeh  (1970)  pointed 
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out  that  Yuen’s  analysis  is  valid  to  second  order,  but  that  his 
third-order  term  is  incorrect  for  wave  numbers  dose  to  the 
cutoff  wave  number.  A  comparison  between  an  experimental 
wave  profile  and  one  calculated  using  Yuen’s  formt.V:on  was 
made  by  Rutland  and  Jameson  (1970);  however,  it  wav  for 
relatively  large  displacements  and  showed  onl>  qualitative 
agreement.  Lafrance  (1975)  used  a  drop  spectrometer  to  meas¬ 
ure  the  size  of  the  main  and  satellite  drops  as  a  function  of 
the  perturbation  frequency.  He  compared  these  to  sizev  cal¬ 
culated  from  a  modified  version  of  Yuen’s  analysis  and  found 
good  agreement  which  is  surprising  since,  in  principle,  the 
theory  is  invalid  in  the  region  where  drops  are  being  formed. 
Chaudhary  (1977)  indicated  that  Lafrance's  analvsis  couldn't 
satisfy  its  own  initial  conditions,  and  Lafrance's  multiple-time 
scale  failed  in  the  third-order  problem.  Consequent!) .  Chaud¬ 
hary  and  Redekopp  (1980a)  and  Chaudhary  and  Maxworthy 
(1980b)  used  a  comparable  model  and  performed  experiments 
on  the  jet  behaviors  and  satellite  drop  formation.  With  a  third- 
order  component  in  the  disturbance,  Chaudhary  et  al.  (19*  \ 
1980a,  1980b)  were  able  to  control  experimentally  the  satellite 
behaviors,  but  his  theory  didn’t  agree  with  his  experiment .  Lee 
(1974)  derived  a  pair  of  nonlinear  one-dimensional  equations 
for  the  jet  radius  and  its  axial  velocity.  He  further  solved  hiv 
one-dimensional  equations  numerically  for  the  infinite  jet 
problem  and  obtained  main  drops  and  satellite  drops.  Boev 
(1978a,  1978b,  1979)  expanded  the  Cosserat  theorv  of  Green 
(1976)  to  develop  a  new  analytical  model  for  higher  harmonics 
based  on  spatial  and  temporal  instability.  The  theory  is  capable 
of  describing  the  satellite  drop  behavior.  The  numerical  in¬ 
vestigation  of  Shokoohi  and  Elrod  (1987)  by  using  the  voniciiy- 
stream  function  method  shows  that  the  nonlinear  terms  are 
responsible  for  the  nonuniformity  of  fhesizeof  drops  produced 
in  jet  disintegration.  The  computations  were  limited  o  a  low 
Reynolds  number  by  numerical  stability  considerations  and  a 
computational  time  requirement.  Phinney  (1972)  carried  out 
experiments  that  Lad  him  to  support  the  conjecture  that  the 
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theory  and  understanding  that  were  developed  in  connection 
with  the  breakup  of  laminar  jets  can  be  used  as  a  guide  for 
turbulent  jets  as  well. 

The  nonlinear  studies  cited  in  the  previous  paragraph  are 
relesant  to  the  ink  jet  breakup.  By  use  of  linear  theory  of 
convective  instability  Lin  and  Lian  (1990)  showed  that  there 
are  two  physically  distinctive  modes  of  jet  breakup.  The  first 
is  the  Rayleigh  mode  w  hich  produces  drops  of  diameters  com¬ 
parable  with  the  jet  diameter.  The  second  is  the  Taylor  mode 
which  atomizes  the  jet  into  droplets  of  diameters  much  smaller 
than  the  jet  diameter.  Lin  and  Creighton  (1990)  showed  that 
the  former  mode  is  caused  by  capillary  pinching,  but  the  latter 
mode  is  caused  by  the  interfacial  pressure  fluctuation. 

To  the  best  knowledge  of  the  present  authors,  this  investi¬ 
gation  is  the  first  nonlinear  analysis  of  the  jet  breakup  in  the 
atomization  regime  taking  into  account  the  viscosities  of  the 
liquid  jet  and  the  ambient  gas  as  well  as  the  gravitational  force. 

Analysis 

We  shall  resort  to  the  application  of  the  nonlinear  theory 
of  hydrodynamic  stability  in  which  disturbances  are  super¬ 
imposed  upon  a  stationary  flow  field.  These  disturbances  are 
time  dependent,  and  it  is  the  aim  of  the  theory  to  look  for 
disturbances  whose  amplitude  increases,  decreases,  or  does  not 
vary  with  time.  Disturbances  of  the  f?rcr  type  we  shall  call 
unstable,  whereas  those  of  the  second  or  third  type  will  be 
called  stable.  We  shall  focus  our  attention  on  disturbances  of 
wave  form  that  repeats  itself  every  wavelength.  Such  an  anal¬ 
ysis  is  called  a  temporal  stability  analysis,  as  opposed  to  a 
spatial  stability  analysis  where  wave  forms  repeat  themselves 
periodically  in  time.  Temporal  stability  analysis  is  adequate 
for  an  infinitely  long  jet  compared  to  the  jet  radius,  while  a 
spatial  stability  analysis  is  adequate  for  a  finite  jet.  Keller  et 
al.  (1973)  have  shown,  in  a  linear  analysis,  that  the  growth 
rate  for  these  two  cases  coincide  when  the  jet  velocity  is  much 
greater  than  the  characteristic  velocity  o/{p )/?).  The  approach 
we  follow  is  similar  to  that  of  Lin  and  Ibrahim  (1990). 

Consider  the  circular  liquid  jet  of  radius  R  that  is  issued 
into  a  viscous  gas  bounded  by  a  cylinderical  wall  at  a  distance 
S  from  the  jet  axis  as  shown  in  Fig.  1 .  Both  the  liquid  and  gas 
are  assumed  to  be  incompressible  and  Newtonian.  The  onset 
of  ligament  formation  and  its  subsequent  breakup  into  droplets 
are  governed  by  the  continuity  and  Navier-Stokes  equations 
of  the  two-dimensional  disturbed  flow  which  can  be  written 
in  dimensionless  form  as 

(  L  +  w)„+  (  U+  u),/r+  (  w)  =0  (/=1,2)  (1) 

(  L'-r  U)„+  (  u+  M),(  L’+  U),r+  (  Vi’),(  U+  u)t: 

=  ~  (P\/pMP  +  p),r  +  *,/(i'iRe)[<l/+  u),„ 

+  ((/^«)(f/r+  (U+  w),::-  {U-u),/r\  (2) 

(  W  4  w),f  +  (  U+  w),(  H  4-  w)(f+ 

=  ~  (pi/p,)(P  +  p),;- Fr~  1  +  (v/p]  Rc)|(H’+ 

+  (W'+  H-)(r/r+  (  W’4-  h  ),..]  (3) 

where  L'and  Hare,  respectively,  th.  radial  and  axial  velocities 
of  the  undisturbed  flow ,  u  and  w  are,  respectively,  the  radial 
and  axial  velocities  of  the  disturbance,  and  P  and  p  are,  re¬ 
spectively,  the  undisturbed  flow  and  the  disturbance  pressures. 
Subscript  1  is  for  liquid  and  2  is  for  gas.  Subscripts  /,  r,  and 
c  stand  for  partial  differentiations  with  respect  to  time,  radial 
distance,  and  axial  distance,  respectively.  Re  is  the  Reynolds 
number  based  on  the  liquid  properties.  Re  =  H'o/?/**, /where 
H'o  is  the  axial  velocity  along  the  jet  axis,  R  is  the  jet  radius, 
and  i'j  is  the  liquid  kinematic  viscosity.  Fr  is  the  Froude  num¬ 
ber,  Fr  =  W  V(gfl),  where  g  is  the  gravitational  acceleration. 
All  lengths  are  made  nondimensional  by  the  jet  radius  R 
Velocities  are  made  dimensionless  b_>  dividing  b>  H0.  Time  is 


Fig.  1  Definition  sketch 


made  dimensionless  by  /?/W'0.  Thv  pressure  is  made  dimen¬ 
sionless  by  dividing  by  piH'5,  where  p\  is  the  liquid  density. 

The  continuity  and  Navier-Stokes  equations  of  the  undis¬ 
turbed  (basic)  flow  are  given  by 

=  0  (4) 

u,  =  0,  P„  =  0  (5) 

0=  Fr’-f  *,/(>•,  Re)(W„,+  »•„//•).  (6) 

Note  that  W  n  =  0  because  the  basic  flow  is  steady.  H  and 
H‘;--  both  vanish  by  virtue  of  the  continuity  equation  which 
also  leads  to  the  conclusion  that  H\  is  a  function  o?  radial 
distance  alone. 

Substitution  from  the  governing  equations  of  the  basic  flow 
into  the  equations,  which  govern  the  disturbed  flow ,  and  elim¬ 
inating  the  pressure  by  cross  differentiation  yield; 

{d  -  i//(i'|Re)£2]£‘V,+  W;£:d,:-rvf!-(  W \/r)r 

=  2E\  i:  -  4 ( r ”  ;E:  v , ),]  C ) 


where  4  is  the  dimensionless  disturbance  stream  function  de¬ 
fined  as, 

^  ,-/rt  w,  =  -  4  Jr  and  E2  =  drr~  ( 1  /r) dr  +  dzz. 

The  solution  of  the  basic  state  equations  subject  to  the  bound¬ 
ary  conditions, 


H'dO)  =  -  1 
H:(i)  =  0 
M  ,(!)=  W\(l) 

H  lr(l)  =  ArH  v(l) 
(P{~P:)  =  \\e 


(finite  velocity  along  the  jet  axis) 

(no-slip  velocit>  at  the  wall) 

(velocity  continuity  at  the  interface) 
(tangential  force  balance  at  the  interface) 
(norma!  force  balance  at  the  interface) 


is  given  by  Ibrahim  (1989)  as 


fL,  =  —  14-  \r /[  V—  ( I  -  /3*)J f  I  —  ((1  —  Q)  (4.\  )] 

x  (Re  Tr)[2  In  8  +  {\  -  tf:)]  1  (8) 

H;=  -  (3:-r)/( A  -  <1~0:)]|  1 

-|(l-0)/(4A)](Re/Fr)|2  In  3  +  (l-d:)]| 

*+  1(1  -  C)/(4.V)l(Re/Fr)l^:-  r  -  2  In  (o  r)]  (9) 


where  3  is  the  aspect  ratio,  3  =  S  7?,  with  S  being  the  distance 
from  the  jet  axis  to  the  wall  and  We  is  the  Weber  number.  We 
=  0/(p,  o  is  the  surface  tension.  .V  is  the  gas  to  liquid 

dynamic  viscosity  ratio  and  Q  is  the  gas  to  liquid  density  ratio. 

The  solution  to  Eq.  (7)  is  sought  in  the  form  of  the  Fourie* 
Bessel  series; 
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Table  1  Demonstration  of  convergence  of  the  eigenvectors. 
Re  =  400,  Re/Fr  =  0.,  We  =  0,0025,  Q  =  0.0013,  K  = 
0.018,  0  =  10.,  k  =  0.3 


l. 

A. 

A*: 

i?(  -  1,  U 

1 

3 

15 

0  001704 

2 

3 

15 

O.OOI818 

3 

3 

15 

0.001 S?4 

2 

2 

8 

0.001797 

tl,(r,z,r)  =  r  F(X,nr)  ^  [a,nl{r)cos(lk:) 

n  *  I  /-0 

+  /?,„/( /)sin(/Ac)]  (10) 
where,  following  Lin  and  Ibrahim  (1990), 

F(\l„r)  =  J]‘  *  ..r)  (11) 

^(X:nr)  =  7,(X>r)y1(X>/S.  -  7,  ( X2n3)  V,  (X2rr)  (12) 

where  Jx  and  are  the  first-order  Bessel  functions  of  the  first 
and  second  kind,  respectively,  k  is  the  wave  number,  and  X<n 
are  obtainable  from  the  orthogonality  conditions  of  F(\inr) 
given  by 

x  i  ( x !  „ )  -  y,  (x)n)  =  o  (13) 

\:nF*  ( X:„ )  -  F{  X:„ )  =  0  (14) 

where  is  the  Bessel  function  of  the  zero  order  of  the  first 
kind,  and 

r  ( x> )  =  J0  ( x> )  y ,  ( x>5 )  -  y,  ( x:„0 )  y0  ( x:r )  ( 1 5) 

where  >'0  is  the  Bessel  function  of  the  zero  order  of  the  second 
kind. 

To  obtain  the  even  mode  of  the  governing  equations  we 
multiply  Eq.  (7)  by  F{\mr)  cos (y‘A*c)  and  integrate  with  respect 
to  r  between  the  limits  e\  and  q>  for  the  liquid  and  the  limits 
e:  and  q:  for  the  gas  where,  e\  =  0,  Q\  =  e:  =  (I  ■+■  17),  q: 
=  3,  and  integrate  with  respect  to  c  between  the  limits  0  and 
k%  where 

riiZJ)  =  \h(ZJ)-\] 

is  the  dimensionless  free  surface  displacem  nt  with  h  being  the 
distance  from  jet  free  surface  to  jet  axis.  Using  a  Taylor  series 
expansion  about  r  =  1  we  can  express  integrations  between 
the  limits  1  and  (1  +  tj)  as, 

A 

«  »  0 

where  \trr  is  anv  function  of  r,  and  A'  is  the  number  of  terms 
in  the  Taylor  series  expansion  required  for  convergence  of  the 
series  within  a  specific  tolerence.  The  solution  of  i?(c,  /)  is 
sought  in  a  Fourier  series  of  the  form 

L 

v (c,/)  =  2  lCtU)co%llkz)+dtU)stn(ikz)].  (16) 

/-o 

The  odd  mode  of  the  governing  equations  is  obtained  by  mul¬ 
tiplying  Eq.  (7)  by  F(X,mr)s\n{jkz)  and  integrating  with  respect 
to  r  between  the  limits  e]  and  q\  for  the  liquid  and  the  limits 
e2  and  q:  for  the  gas  and  integrating  with  respect  to  z  between 
the  limits  0  and  2r/k.  Here,  j  »  0,  1,  2,  . . . ,  L  and  m  =  1, 
2,  . . . ,  A',  for  liquid  and  m  =  1,  2,  ...»  As  for  gas,  where 
A’,,  A2,  and  L  are  the  number  of  terms  in  Eq  (10)  required 
to  yield  a  sufficiently  accurate  representation  of  $,{r9zj).  The 
reason  behind  including  only  (A’:  -  5)  equations  of  the  gov¬ 
erning  equations  of  the  gas  is  that  only  five  of  the  eight  bound- 
ar>  conditions  arc  not  readily  satisfied  by  the  series 
representation  in  Eq.  (10)  and  that  we  are  adopting  the  Lanczos 
(1956)  tau  method  in  our  solution  by  the  Spectral  method.  The 
three  boundary  conditions  that  are  readily  satisfied  by  Eq  (10) 
are  those  which  require  finite  radial  and  axial  components  of 
disturbance  velocity  along  the  jet  axis  and  the  vanishing  of  the 


radial  component  of  velocity  at  the  wall.  The  remaining  five 
boundary  conditions  are  those  w  hich  require  zero  axial  \  elocity 
at  the  wall,  continuous  radial  and  axial  velocity,  tangential 
stresses  and  normal  stresses  at  the  liquid-gas  interface,  which 
can  be  written  in  equation  form  as, 

ws(/3.r)  =  0  (17) 

u](h>z)  =  u2(h,z)  (18) 

W\{h)  +  wj(/j)  =  -»■  h:( A)  (19; 

I u]:(h,z)  +  W  lr ( * )  +  mv(A,c)] 

=  A] uZz(h%z)  +  W ’2,(h)  t  h :,(A,c)]  (20) 
-  l-[Pi(h)+Pi(htz)}  +  <2/Re)i/„(Atc)| 

+  {  -  lP2(h)  +p2(h,z)]  -r  (2A7Re)i/;,(/i,c)  I 

-We(l-i)  +  jj‘-i);:-rj:/2)  =  0,  (21) 

respectively.  Note  that  the  last  term  in  Eq.  (21).  which  rep¬ 
resents  the  surface  tension  stresses,  is  simplified  to  its  second- 
order  form  in  disturbance  components  (Ibrahim,  1989)  We 
are  going  to  adopt  only  second-order  analysis  in  the  present 
work  due  to  the  complexity  of  the  equations.  The  boundary 
conditions  in  Eq.  (18)  through  Eq.  (21)  apply  at  the  liquid-gav 
interface  which  is  not  Know  n  a  priori.  These  boundary  con¬ 
ditions  are  expanded  in  the  Taylor  series  of  second  order. 
Carrying  out  this  expansion  and  using  the  basic  state  boundary 
conditions,  Eq.  (17)  through  Eq.  (21)  become,  in  terms  of  the 
disturbance  stream  function, 


*!;,<0.O  =  O 

(22) 

-  [V2r.-(1.C)  -  1 

.0)1 

(23) 

z)))[ 

(24) 

(0  M,H‘ 

T-;  ~  ->■  vi'i-n;-  C:;)  -  »)([  -(Z^v  ),). 

+  (v.r-  Cr))  )];  =  (!  -  Q)  (Re  Frt<n 

(25) 

[(P,  P;)Re|(’Jr;-  TJV-.,-  T)X:;)  -  M  ,(  1  )  V  -  —  V;(v,  -  V, 

-) 

+  (  H  (  1  )  -  V  r)  V  -  W  (  1  >77  V  -  Wd  )7):V 

1 

+  (P.  •>  i )  1  (  -  (D: V  )  r]r  ,  1  -  (  3 V  r-:  -  2C  ::  ) 

::  V  rr--  —  4v  ::)  V;”;]  "*■  r)  -  ( {  3  V  ~  V  rr;  ”■  4v 

;)**  V  • 

J!): 

Re  V\  e( 77 _  +  7]...  -  21777.  -  rj: 

1  =  0, 

(26) 

where 

is  the  dynamic  viscosity.  An  additional 

equat: 

ion  iv 

needed  to  account  for  the  unknown  r}.  This  equation 

is  obtained 

from  the 

relation  between  the  radial  component  of 

veloc 

ity  of 

the  liquid  at  the  liquid-gas  interface  u  (h%z)  and  i?; 

u]  ( h,z )  =  D-n  Dt. 

(2**) 

The  boundary  conditions  corresponding  to  the  even  and  odd 
modes  are  obtained  by  multiplying  Eq.  (22)  through  Eq  (2“> 
by  cosf/A'c)  and  $in(jkz),  respectively,  and  integrating  with 
respect  to  z  between  the  limits  0  and  2 z/k. 

The  number  of  unknowns  is  equal  to  the  number  of  equa¬ 
tions  is  equal  to  2(1  +  l)(Aj  +  A:)  +  2(L  +  1 ).  The  problem 
is  solved  as  an  initial-value  problem  with  the  initial  condition'' 
given  by  the  solution  of  the  corresponding  linear  problem  av 
obtained  by  Lin  and  Ibrahim  (1990).  A  fourth-order  mulustep 
predictor  corrector  scheme  described  in  Burden  and  Faires 
(1985)  is  used.  The  values  of  A’j,  A:,  and  L  are  chosen  to  be 
3,  15,  and  2,  respectively,  as  justified  by  the  results  of  the 
convergence  test  shown  in  Table  1 .  More  details  on  the  solution 
procedure  can  be  found  in  Ibrahim  (1989). 

Results  and  Discussion 

Figure  2  depicts  a  comparison  between  the  linear  and  non¬ 
linear  predictions  of  i  mporal  evolution  of  the  free  surface 
displacement  in  the  atomization  mode.  The  linear  solution  iv 
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Fig  2  Comparison  of  the  linear  and  nonlinear  predictions  of  the  tem 
poral  evolution  of  the  interface  displacement  in  the  atomization  mode 


Fig  3  Comparison  of  the  linear  and  nonlinear  predictions  of  the  in¬ 
terfacial  profile  in  the  atomization  mode 


being  that  of  Lin  and  Ibrahim  (1990).  It  is  seen  that  the  non¬ 
linear  analysis  predicts  a  peak  that  occurs  sooner  and  reaches 
a  higher  magnitude  than  that  of  the  linear  analysis.  In  other 
\*ords.  the  initial  temporal  evolution  of  finite  disturbances 
occurs  at  a  higher  growth  rate  and  frequencs  than  those  as¬ 


sociated  with  infinitesimal  disturbances.  This  means  that  finite 
disturbances  are  more  unstable  than  infinitesimal  disturbances 
which  leads  to  the  conclusion  that  if  linear  theory  predicts  that 
a  jet  is  unstable  to  infinitesimally  small  disturbances,  then 
nonlinear  theorv  further  predicts  that  the  disturbances  will 
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Fig  4  Nonlinear  interfacial  profile  in  the  capillary  mode.  f(-  1,  0)  = 
0.00146 


continue  to  grow  after  they  become  finite  and  ultimately  the 
jet  will  be  atomized. 

Figure  3  describes  the  interfacial  profile  of  the  jet  surface 
in  the  atomization  mode  as  predicted  by  both  the  linear  and 
nonlinear  analyses.  The  horizontal  axis  is  normalized  with 
( -  k  f2x)  so  that  the  interfacial  profile  is  depicted  for  one 
wavelength  2 rJk.  Since  the  present  analysis  is  temporal,  this 
profile  repeats  itself  every  wavelength  2ir /k.  The  negative  sign 
in  the  normalization  (-k/ 2*)  is  to  make  the  horizontal  co¬ 
ordinate  increase  in  the  flow  direction  since  the  axial  coordinate 
Z  is  chosen  to  be  negative  downward  as  shown  in  Fig.  1 .  Figure 
3  shows  that  the  nonlinear  analysis  produces  larger  defor¬ 
mation  of  the  jet  surface  than  the  linear  instability  analysis. 
This  would  eventually  result  in  larger  droplet  size,  after  dis¬ 
integration  of  the  jet,  than  those  that  would  be  predicted  by 
the  linear  analysis. 

Figure  4  shows  the  configuration  of  the  liquid-gas  interface 
in  the  capillary  mode  of  instability  as  predicted  by  the  nonlinear 
analysis.  It  is  seen  that  the  initially  sinusoidal  wave  develops 
into  an  asymmetric  wave  with  a  higher  harmonic.  This  higher 
harmonic  was  discovered  in  the  analysis  of  Yuen  (1968)  and 
the  experiments  conducted  by  Rutland  and  Jameson  (1970), 
Lafrance  (1975),  and  Taub  (1976)  among  others.  The  asym¬ 
metry  is  due  to  the  higher  harmonics  as  well  as  feedback  into 
the  fundamental.  Since  the  present  analysis  is  only  of  second 
order,  we  don’t  expect  to  see  more  than  one  swelling  between 
the  crests  of  the  primary  disturbance  waves.  As  a  matter  of 
fact,  this  is  exactly  what  Rutland  and  Jameson  (1970)  and 
Lafrance  (1975)  observed  up  to  the  point  of  jet  breakup.  After 
breakup,  of  course,  the  ligaments  between  large  drops  can 
often  disintegrate  further,  in  which  case  many  undulations  per 
primary  wavelength  could  be  observed  as  reported  by  Rutland 
and  Jameson  (1970).  However,  Taub  (1976)  reported  up  to 
three  undulations  before  jet  breakup.  It  should  be  noted  that, 
our  results  are  obtained  at  a  much  smaller  time  than  those  of 
other  authors  and  hence  the  amplitude  of  the  interface  dis¬ 
placement  in  Fig.  4  is  much  smaller.  Unfortunately,  the  CPU 


time  required  to  advance  the  solution  one  dimensionless  time 
unit  was  too  large  to  attempt  to  obtain  results  at  a  larger  time. 

The  occurence  of  higher  harmonics  explains  the  common 
observation  that  once  the  jet  breakup  it  disintegrates  inio  drop¬ 
lets  of  various  sizes.  This  is  in  contrast  to  the  predictions  of 
the  linear  theory  of  jet  disintegration  into  droplets  of  uniform 
size. 

Referring  back  to  Fig.  3,  it  is  obvious  that,  unlike  Rayleigh 
mode,  the  results  don’t  indicaie  any  higher  harmonics  for  the 
atomization  mode.  This  could  be  explained  by  observing  that 
the  amplitude  of  the  interface  displacement  in  the  capillary 
mode  is  two  order  of  magnitudes  larger  than  in  the  atomization 
mode  so  it  is  easier  for  internodal  undulation  of  finite  ampli¬ 
tude  to  take  place  in  the  capillary  mode.  Also,  the  higher 
growth  rates  and  frequencies  associated  with  the  atomization 
mode  could  be  the  reason  why  it  is  improbable  for  asymmetry 
to  set  in.  However,  it  is  possible  that  if  the  nonlinear  solution 
is  carried  out  for  a  larger  time,  we  will  see  higher  harmonics 
for  the  atomization  mode.  Such  an  investigation  was  not  car¬ 
ried  out  because  of  the  enormous  CPU  time  required. 


Conclusions 

While  no  nonlinear  stability  was  found  for  unstable  linear 
disturbances,  stable  linear  disturbances  were  found  to  be  also 
stable  with  respect  to  weakly  nonlinear  disturbances  in  a  wide 
range  of  flow  parameters.  The  nonlinear  instability  of  the 
atomization  mode  indicated  that  finite  disturbances  are  more 
unstable  than  infinitesimally  small  disturbances.  No  higher 
harmonics  are  observed  for  the  atomization  mode  for  the  time 
range  investigated. 

Asymmetry  of  the  interfacial  profile  is  revealed  in  the  non¬ 
linear  analysis  of  capillary  instability.  This  explains  the  ob¬ 
served  satellite  formation  in  ink  jets.  The  amplitude  of  the 
interfacia]  disturbances  are  two  order  of  magnitude  larger  than 
those  associated  with  the  atomization  mode  of  instability.  Lin- 
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ear  theories  are  adequate  for  predicting  the  onset  of  either  the 
Rayleigh  or  the  Taylor  mode  instability.  The  outcome  of  atom¬ 
ization  process  requires  fully  nonlinear  treatment. 
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The  effect  of  compressibility  of  fluids  on  the  atomization  of  a  liquid  jet  emanating  from  a  nozzle  into  the 
ambient  gas  is  investigated  by  use  of  linear  stability  analysis.  An  increase  in  the  compressibility  of  the  liquid 
gas  is  shown  to  promote  the  generation  of  smaller  droplets,  although  at  a  slower  rate.  It  is  shown  that  in  some 
parameter  ranges,  a  small  change  in  compressibility  may  result  in  a  qualitative  change  of  the  mechanism  of 
atomization. 


Introduction 

TOM1ZATION  is  a  process  of  breaking  up  a  liquid  jet 
into  droplets  of  diameter  much  smaller  than  the  jet  di¬ 
ameter.  This  process  is  widely  used  in  industrial  applications 
including  fuel  injections  in  internal  combustion  engines.  A 
good  knowledge  of  the  fundamental  mechanism  of  atomiza¬ 
tion  is  essential  for  raising  the  combustion  efficiency  and  re¬ 
ducing  the  environmental  pollution.  Our  understanding  of 
atomization  is  not  yet  complete. 

A  concise  review  of  experimental  and  theoretical  work  on 
atomization  was  given  by  Reitz  and  Bracco.12  A  more  ex¬ 
tensive  exposition  of  the  subject  can  be  found  in  the  work  of 
Lefebvre/  The  atomization  process  as  the  consequence  of  the 
aerodynamic  instability  of  a  jet  emanated  into  ambient  gas 
was  investigated  by  Taylor/  Levich/  and  others.  They  con¬ 
sidered  instability  with  respect  to  temporally  growing  dis¬ 
turbances.  Atomization  as  the  consequence  of  convective  in¬ 
stability  caused  by  spatially  growing  disturbances  was  studied 
by  Lin  and  Kang/  Lin  and  Lian.7  and  Lin  and  Creighton/ 
In  all  of  the  theoretical  works  mentioned  above,  the  fluids 
are  assumed  to  be  incompressible.  On  the  other  hand  in  many 
applications,  the  jet  speed  is  such  that  the  Mach  number  of 
the  ambient  gas  relative  to  the  jet  may  reach  transonic.  For 
future  applications  it  may  even  reach  supersonic.  Therefore, 
the  compressibility  of  the  fluids  must  be  taken  into  account. 

Here,  we  investigate  the  effects  of  compressibility  on  atom¬ 
ization.  The  problem  of  instability  of  an  inviscid  compressible 
jet  emanating  into  an  inviscid  compressible  gas  is  formulated  in 
the  next  section.  The  disturbance  is  allowed  to  grow  both  tem¬ 
porally  and  spatially .  The  results  of  the  solution  to  the  formu¬ 
lated  problem  are  used  to  demonstrate  that  the  compressibility 
brings  about  not  only  quantitative  but  also  qualitative  differences 
in  the  mechanism  of  atomization  at  its  onset. 

Formulation  of  Problem 

Consider  the  onset  of  instability  of  a  jet  of  inviscid  com¬ 
pressible  fluid  emanating  from  a  circular  cylindrical  nozzle  of 
radius  a  into  an  unbounded  inviscid  compressible  fluid.  The 
governing  equations  of  motions  are 

pXb.  +  V,T)V,  =  -XP,  (1) 
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to,  +  V-(p y,)  =  0  (1  =  1.2)  (2) 

where  the  subscripts  /  =  1.2  denote  respectively  the  jet  fluid 
and  the  ambient  field,  p  is  density,  t  is  time.  V  is  the  velocity 
vector,  and  P  is  the  dynamic  pressure.  The  basic  state  that 
satisfies  Eqs.  (1)  and  (2)  exactly  in  the  absence  of  gravity  is 
given  by 

V,  =  (0,.  1/„  U\)  =  (0,  0,  We).  p,  =  constant 

P,  =  constant.  (0  s  R  s  a) 

=  (&;,  ^j.  ft.)  =  (0.  0.  0).  p,  =  constant 

P:  =  Px  -  S/a.  (a<  R  <*) 

where  the  overbar  denotes  the  basic  state  quantities.  U ,  V, 
and  IV  are  respectively .  the  radial,  azimuthal  and  axial  com¬ 
ponents  of  velocity  in  the  cylindrical  coordinate  system  {R. 
0 .  Z).  Wtt  is  the  constant  jet  speed,  and  5  is  the  surface  tension. 
Z  is  positive  in  the  flow  direction.  The  basic  state  liquid  jet 
is  of  constant  radius  with  smooth  surface.  Thus  neither  shock 
waves  nor  Mach  waves  exist  in  the  gas  and  the  liquid  jet,  and 
the  densities  in  both  fluids  are  constant  in  the  basic  state. 

We  investigate  the  stability  of  this  basic  state  with  respect 
to  the  disturbance  defined  by 

v,  =  v,  +  v;.  p,  =  p,  +  p p,  =  p,  +  p;  (3) 

where  primes  are  used  to  designate  the  perturbations  of  the 

variables.  Substituting  Eq.  (3)  into  Eqs.  ( 1 )  and  (2).  retaining 
only  linear  terms,  we  have 

(b,  +  6„W..*/)V;  =  -VP,  (4) 

(8,  +  6„w;,d^)p;  =  -p(v  17  (5) 

where  is  1  or  0  depending  on  if  i  =  1  or  i  *  1.  Taking 
the  divergence  of  Eq.  (4)  and  then  using  Eq  (5).  we  have 

<*,  +  6„vv  =  v-p;  <6) 

P]  and  p\  are  related  locally  by  the  speed  of  sound  c;  = 

dP'Jdp Thus  Eq.  (6)  can  be  written  as 

(A,  +  6uw„a,vp;  =  rfv-p;  (7) 

The  corresponding  linearized  boundary  conditions  for  Eq  (7) 
at  the  interface  R  ~  a  +  d  are  the  dynamic  force  balance 

PJ  -  P;  =  S{a  :  -  a  -  b,,)d  (S) 
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and  the  kinematic  conditions 

(a,  +  6uw„yd  « 


For  the  special  case  of  an  incompressible  jet  in  a  vacuum. 
Q  =  0 ,  c,  — *  — *  0,  thus  A,  — *  A,  A:  — »  sc  and  Eq.  (16) 

(9)  is  reduced  to,  if  w  remains  finite 


where  d  is  the  displacement  of  the  interface  from  the  basic 
state  cylindrical  surface.  Note  that  P\  in  Eqs.  (8)  and  (9)  are 
to  be  applied  at  R  =  a. 

The  normal  mode  solution  of  the  differential  system  Eqs. 
(7),  (8),  and  (9)  will  be  sought  in  the  form 


X 

P,W;,pM 

p: 

= 

P,Wip:(r) 

d 

of  J 

•exp[/(Az  +  nS  -  o>r)] 


(10) 


where  (r,  z,  /).  and  r  are  dimensionless  variables  defined  by 
(r,  z,/)  =  (fl.  Z,  d)a~K  and  r  =  t(WJa);  and  (A,  n)  and  u> 
are  the  dimensionless  wavenumber  vector  and  wave-fre¬ 
quency,  respectively. 

Substitution  of  Eq.  (10)  into  Eqs.  (7),  (8).  and  (9)  yields 
the  governing  equation  for  the  dimensionless  perturbation 


pressure  amplitudes 

{£):  +  r~xD  -  \rt:r~:  +  A:  -  M](k  -  <o):]}p,  =  0  (11) 

{D:  +  r  'D  -  (n->-  +  A-  -  M:u>2]}p2  =  0  (12) 

and  the  corresponding  boundary  conditions  at  r  =  1 

Pz  -  Pi  +  W,(A2  +  -  1)/  *  0  (13) 

(A  -  w)~f  -  p\  =  0  (14) 

Q(o-f-p’2  =  0  (15) 

where 

VV  =  5/p,W;,a,  Q  =  p:/p,,  Af,  =  WJct. 

M  =  WJc2 ,  D  =  a/dr 


Note  that  VV,  is  the  Weber  number,  M,  is  the  Mach  number 
of  the  jet,  and  M  is  the  Mach  number  of  the  ambient  fluid 
with  respect  to  a  frame  of  reference  moving  with  the  jet  at  a 
speed  W„.  The  bounded  solution  of  Eqs.  (11)  and  (12)  are 

p,  =  v4/f#(Atr),  A,  =  [A2  -  M-(k  -  w)2}': 

p:  =  BK„(A:r),  A:  =  [A2  -  MV),: 

where  A  and  B  are  the  integration  constants  and  /„  and  K„ 
are  respectively  the  modified  Bessel  function  of  the  first  and 
second  kinds.  Substitution  of  these  solutions  into  Eqs.  (14) 
and  (15)  yields 


A  =  (A  -  «u)-///:(A1)A, 

B  =  Qu2f!KjA2)  A, 

where  primes  now  denote  the  differentiation  <f/</(A,r).  It  fol¬ 
lows  from  Eq.  (13)  that 

=  (A  -  <u)-7„(A,)  _  <?orX(A:) 

"  A,/;, (A,)  a:/c;,(a:) 

-  VV(A:  +  n-  -  1)  =  0  (16) 

For  a  given  set  of  four  flow  parameters  VV  .  Q .  Af,.  and  M. 
this  is  the  secular  equation  the  solution  of  which  gives  the 
characteristic  wave  frequency  <a  and  the  wavenumber  A. 


IJk)(,k{k)  —  +  *()  -  =  o  (i7) 

This  is  the  characteristic  equation  used  by  Keller  et  al M  and 
Leib  and  Golstein10  in  their  studies  of  the  spatial  instability 
of  the  Rayleigh  jet.  For  the  special  case  of  an  incompressible 
jet  emanating  into  an  incompressible  fluid,  M .  M,  — ►  0.  thus 
Aj  =  A:  =  A,  and  Eq.  (16)  is  reduced  to,  if  w  remains  finite: 


(A  -  <oy/n(k) 
K(k) 


Q^Kn(k) 


-  A(A2  +  rt2  -  1)VV,  =  0 
(18) 


This  equation  agrees  with  Eq.  (1)  of  Lin  and  Lian7  when  the 
Reynolds  number  of  their  equation  is  let  to  approach  infinity 
and  their  o>  and  A  are  respectively  replaced  by  -io>  and  -  A. 
This  change  of  variable  is  necessary  because  of  the  different 
positive  axial  directions  and  the  different  definitions  of  the 
wave-frequency  which  are  used  in  these  two  papers.  When 
the  surface  tension  is  put  to  zero,  Eq.  (16)  reduces  to  the 
characteristic  equation  corresponding  to  that  of  Tam  et  al.n 
except  their  equation  is  in  terms  of  the  Hankel  function  and 
the  Bessel  function  of  the  first  kind.  It  is  known  that  the 
axisymmetric  disturbance  grows  the  fastest  in  an  incompres¬ 
sible  jet  in  the  absence4  or  presence12  of  the  ambient  field. 
Here  we  assume  that  the  same  is  true  for  a  compressible  jet. 
Hence  we  treat  the  case  of  n  =  0  in  Eq.  (16).  Here  both  w 
=  <or  +  itj,  and  A  =  A,  4  ik,  are  complex.  Their  imaginary 
parts  w ,  and  A ,  give  respectively  the  temporal  and  the  spatial 
growth  rates.  Their  real  parts  are  the  wave  frequency  and  the 
wave  number  respectively. 


Results 

In  order  to  demonstrate  the  significant  effect  of  compress¬ 
ibility  on  the  atomization  of  liquid  jets,  we  present  first  the 
results  obtained  with  the  assumption  that  the  fluids  involved 
are  incompressible.  Figure  1  shows  the  dependence  of  the 
spatial  amplification  rate  A,  on  the  wave  number  k,  for  the 
flow  parameters  shown  in  the  figure  caption.  A/,  =  M  =  0 
implies  that  both  the  liquid  jet  and  the  ambient  gas  are  in¬ 
compressible.  The  values  of  VV,  =  10~5  and  Q  =  0.0013 
correspond  to  the  case  of  a  water  jet  emanated  at  a  speed  of 
201  m's  from  a  nozzle  of  radius  0.18  mm  into  room  temper¬ 
ature  air  under  1  atmospheric  pressure.  These  values  are  en- 


Fig.  1  Convective  instability  of  incompressible  jet  Q  =  0.0013,  H, 
=  0.00001,  M,  =  0.0,  M  =  0.0. 
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countered  in  many  common  applications.  Along  the  curve  in 
this  figure,  w,  «  0.  It  is  seen  from  Eq.  (10)  that  the  negative 
value  of  k,  signifies  spatial  growth  of  the  disturbance  ampli¬ 
tude  in  the  flow  direction,  and  a>,  =  0  signifies  zero  temporal 
growth  of  the  disturbance.  Thus  the  jet  under  the  specified 
parameter  condition  exhibits  convective  instability.13 14  How¬ 
ever,  these  results  were  obtained  with  the  assumption  that 
both  the  liquid  jet  and  the  ambient  gas  are  incompressible. 
Actually  both  fluids  are  compressible,  and  the  Mach  number 
of  the  above  water  jet  flowing  at  a  speed  of  201  m/s  is  0.135. 
The  corresponding  convective  Mach  number  of  the  ambient 
air,  i.e.,  the  Mach  number  of  the  air  relative  to  a  reference 
frame  moving  with  the  jet,  is  0.59.  With  these  values  of  Mach 
numbers  and  the  rest  of  the  flow  parameters  remaining  the 
same  as  that  in  Fig.  1,  we  obtained  the  solution  curves  of  the 
dispersion  equation  (16).  The  results  are  presented  in  Fig.  2. 
There  are  three  branches  of  the  spatial  growth  curve  =  0. 
The  constant  frequency  curves  wt  =  2.845  and  wr  =  2.866 
being  initiated  from  the  upper  half  plane  approach  the  left 
and  right  branches  respectively  as  the  values  of  w,  are  de¬ 
creased  from  some  positive  values  toward  zero.  Along  the 
left  and  right  branches,  wr  decreases  and  increases  respectively 
with  k The  group  velocities  of  disturbances,  being  given  by 
the  derivation  of  a>,  with  respect  to  k„  are  in  the  downstream 
and  upstream  directions  respectively  for  the  left  and  right 
branches.  Note  the  existence  of  a  saddle  point  with  the  value 
of  u>,  lying  between  2.845  and  2.866.  At  this  point  D0  =  dD0 / 
dk  =  0,  and  the  group  velocity  is  zero.  According  to  the 
theories  of  Briggs13  and  Bers,14  the  appearance  of  the  saddle 
point  signifies  that  the  integration  path  in  the  Fourier  integral 
representation  cannot  be  analytically  deformed  continuously 
until  to,  =  0  is  reached  from  above.  The  disturbance  has  to 
grow  simultaneously  in  time  and  space,  according  to  the  the¬ 
ory  of  Briggs  and  Bers.  Physically,  while  the  atomization  starts 
at  some  distance  downstream  of  the  nozzle  in  convective  in¬ 
stability,  the  jet  breakup  start  right  from  the  nozzle  in  absolute 
instability.  A  more  precise  quantitative  physical  consequence 
of  absolute  instability  requires  the  solution  of  the  initial  value 
problem  of  nonlinear  instability.  A  comparison  of  Fig.  1  with 
Fig.  2  clearly  demonstrates  the  importance  of  the  compress¬ 
ibility  effect  on  the  onset  of  atomization.  It  brings  about  a 
qualitative  difference  in  the  mechanism  of  atomization. 

Absolute  instability  can  be  eliminated  by  increasing  the 
compressibility.  When  the  Mach  numbers  in  Fig.  2  are  quad¬ 
rupled  while  the  rest  of  parameters  remain  fixed,  the  saddle 
point  is  removed  and  the  solution  of  Eq.  (16)  yields  the  am¬ 
plification  curve  of  convectively  unstable  disturbances  shown 


Fig.  2  Absolute  instability  of  subsonic  Jet  Q  *  0.0013,  H  r  *  0.00001, 
M%  *  0.135,  M  *  0.59. 


Fig.  3  Convective  instabibrv  at  subsonic  jet  Q  =  0.0013,  W  ,  =  0.00001, 
M ,  =  0.52,  M  -  2.3 


in  Fig.  3.  Note  that  k,  becomes  zero  beyond  the  critical  wave 
number  kr  -  34,000.  Thus  the  smallest  atomized  droplet  ac¬ 
cording  to  the  linear  theory  will  have  a  radius  of  approxi¬ 
mately  9  ftm.  This  estimate  is  obtained  by  using  the  assump¬ 
tion  that  a  droplet  is  formed  by  capillary  pinching  at  two 
consecutive  troughs  of  the  liquid  jet.  This  estimate  therefore 
tends  to  over  predict  the  diameter  since  the  droplets  are  ac¬ 
tually  formed  from  small  ligaments  branching  out  from  the 
jet.  The  conservation  of  volume,  then  requires  the  following 
relation  between  the  radius  of  droplet  R  and  the  wavelength 
of  disturbance  A: 

tto’A  =  4ttK3/3  (19) 

where  A  is  related  to  the  wave  number  kr  by  k,  =  liraik.  The 
estimated  smallest  droplet  can  be  compared  with  that  ob¬ 
tained  by  neglecting  the  compressibility  effect.  Referring  back 
to  Fig.  1 ,  we  see  that  the  critical  wave  number  above  which 
the  value  of  k,  becomes  zero  is  approximately  130.  By  use  of 
the  same  method,  the  smallest  atomized  droplets  are  esti¬ 
mated  to  be  approximately  60  Thus  the  smallest  droplet 
predicted  by  the  theory  without  taking  into  account  the  com¬ 
pressibility  effect  is  almost  one  order  of  magnitude  larger  than 
that  when  the  compressibility  effect  is  duelv  taken  into  ac¬ 
count.  If  one  reasonably  assumes  that  the  population  of  the 
droplets  of  a  given  diameter  in  a  spray  is  proportional  to  the 
amplification  rate  -kn  then  a  simple  comparison  of  Figs.  1 
and  3  clearly  indicates  that  the  size  distributions  in  the  sprays 
of  compressible  and  incompressible  fluids  are  quite  different. 
An  incompressible  jet  tends  to  produce  a  smaller  range  of 
droplet  diameters  but  with  a  flatter  distribution.  It  should  be 
pointed  out  that  although  a  high  speed  jet  with  M ,  =  0.52  is 
seldom  encountered  in  present  day  practice,  it  will  probably 
be  encountered  in  the  future  in  high  speed  atomization  for 
small  droplets. 

In  the  remainder  of  this  work,  we  present  the  results  dem¬ 
onstrating  the  effect  of  surface  tension,  ambient  gas  density, 
and  jet  speed  on  atomization.  The  values  of  flow  parameters 
for  each  amplification  curve  are  given  in  the  figure  caption. 
When  the  surface  tension  of  the  jet  whose  disturbance  am¬ 
plification  curve  is  given  by  Fig.  3  is  increased  by  a  factor  of 
10,  the  amplification  rate  is  reduced  by  a  factor  of  approxi¬ 
mately  ten  as  shown  in  Fig  4.  This  indicates  that  surface 
tension  is  stabilizing.  The  critical  wave  number  is  also  reduced 
by  a  factor  of  approximately  10,  and  thus  the  diameter  of  the 
smallest  size  of  droplet  obtainable  is  increased  by  a  factor  of 
10* 3  (c.f.  Eq.  (19)).  When  the  Q  in  Fig.  3  is  increased  by  a 
factor  of  10,  the  amplification  rate  is  increased  by  a  factor  of 
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Fig.  4  Surface  tension  is  stabilizing  Q  =  0.0013,  H,  *  0.0001,  Af, 
*  0.52,  M  *  2.3. 


Fig.  5  Ambient  gas  promotes  atomization  Q  =  0.013,  W  ,  =  0.00001, 
A# ,  *  0.52,  M  =  2.3. 


approximately  4  as  shown  in  Fig.  5.  The  critical  wave  number 
is  only  slightly  increased.  However,  the  amplification  curve 
is  considerably  flattened.  Thus  the  effect  of  increasing  the 
ambient  gas  density,  by  increasing  the  gas  pressure  for  ex¬ 
ample.  is  to  enhance  the  atomization  rate  at  the  expense  of 
getting  a  broader  size  distribution.  When  the  Mach  number 
of  the  jet  for  Fig.  3  is  increased  from  0.52  to  a  supersonic 
value  of  1.25  with  M  fixed  at  2.3,  the  amplification  rate  is 
slightly  reduced  and  the  critical  wave  number  is  slightly  in¬ 
creased  as  shown  in  Fig.  6.  For  a  given  gas,  this  signifies  that 
if  the  jet  speed  is  kept  constant  but  the  speed  of  sound  is 
reduced  by  a  factor  of  0.51/1.25,  one  will  observe  a  slight 
reduction  in  the  rate  of  atomization.  Thus,  an  increase  in 
compressibility  of  the  jet  softens  the  atomization  process,  and 
the  most  populous  droplet  tends  to  have  a  slightly  smaller 
diameter. 

Discussion 

The  Reynolds  number  corresponding  to  the  examples  given 
in  Figs.  2  and  6  ranges  from  20,000  to  75,000  if  the  jet  liquid 
is  water.  This  justifies  the  neglect  of  viscosity  in  the  first 
approximation.  Only  the  results  for  the  axisymmetric  mode 
are  presented  here.  This  is  because  the  axisymmetric  mode 


Fig.  6  Compressibility  of  jet  reduces  atomization  rate  Q  =  0.0013, 
H,  *  0.00001,  Af,  =  1.25,  M  =  2.3. 


is  the  most  unstable  mode.  This  can  be  seen  from  Eq.  (16). 
The  first  two  terms  in  this  equation  arise  from  the  pressure 
fluctuation  which  is  destabilizing.  The  W,  term  arises  from 
the  capillary  effect.  The  first  two  terms  in  the  parenthesis  are 
positive  and  represent  respectively  the  effect  of  stabilization 
by  surface  tension  associated  with  the  curvatures  along  the 
axial  and  azimuthal  directions  respectively.  The  only  negative 
term  in  the  parenthesis  arises  from  the  destabilizing  capillary 
pinching  perpendicular  to  the  axial  direction  Thus  the  non- 
axisymmetric  modes  associated  with  n  ^  0  are  more  stable 
The  thin  boundary-layer  effect  near  the  liquid-gas  interface 
was  discussed  by  Lin  and  Lian7  and  Lin  and  Ibrahim15  for  the 
case  of  incompressible  fluids.  The  viscosity  was  found  to  be 
responsible  for  generating  larger  drops  which  scale  with  shear 
waves.  The  wavelength  of  shear  waves  is  several  orders  of 
magnitude  larger  than  the  gas-capillary  length  5/p:U;  = 
W,a/Q.  Note  that  the  wave  length  of  the  most  amplified 
disturbance  in  Fig.  1  scales  with  this  capillary  length,  since 
the  wavelength  A  =  2n alkr.  On  the  other  hand,  the  wave¬ 
lengths  of  the  most  amplified  disturbances  given  in  Figs.  2  to 
6  for  compressible  jet  scale  with  the  liquid-capillary  length 
S/p,W;  =  W^,  this  can  be  verified  easily  from  the  result 
given  in  these  figures.  However,  the  length  scale  of  the  shear 
waves  mentioned  above  was  obtained  with  a  theory  in  which 
the  compressibility  effect  is  totally  neglected.  A  theory  which 
takes  into  account  both  the  effect  of  viscosity  and  compress¬ 
ibility  is  not  yet  available.  A  nonlinear  initial  value  problem 
for  the  elucidation  of  absolute  instability  in  an  atomizing  jet 
is  yet  to  be  formulated. 
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The  absolute  and  convective  instability  of  a  jet  of  an  inviscid  compressible  fluid  emanating  into 
another  inviscid  compressible  fluid  is  investigated.  It  is  shown  that  the  compressibility  of  the 
ambient  fluid,  as  well  as  the  surface  tension,  tends  to  promote  absolute  instability.  The  absolute 
instability  can  be  reduced  to  the  convective  one  by  reducing  the  surface  tension  and/or 
increasing  the  compressibility  of  the  jet  fluid,  depending  on  the  flow  parameters. 


L  INTRODUCTION 

The  capillary  instability  of  an  infinitely  long  liquid  jet 
with  respect  to  temporally  growing  disturbances  was  ana¬ 
lyzed  by  Rayleigh. 1  Keller  et  al?  examined  the  capillary  in¬ 
stability  of  a  semi-infinite  liquid  jet  with  respect  to  spatially 
growing  disturbances.  They  found  that  the  temporal  and 
spatial  disturbances  are  analytically  related  if  the  Weber 
number  is  sufficiently  small.  When  the  Weber  number  is 
sufficiently  large,  they  also  found  faster  growing  spatially 
amplifying  disturbances  that  are  unrelated  to  the  Rayleigh 
mode  of  disturbances.  Leib  and  Goldstein3  showed  that  the 
new  mode  found  by  Keller  et  al  actually  cannot  be  reached 
by  a  given  initially  infinitesimal  disturbance  in  the  sense  of 
Briggs4  and  Bers.5  Using  the  dispersion  equation  obtained  by 
Chandrasekhar,6  Leib  and  Goldstein7  obtained  the  critical 
Weber  number  above  which  an  incompressible  viscous  jet  is 
absolutely  unstable  and  below  which  it  is  convectively  unsta¬ 
ble  as  a  function  of  the  Reynolds  number.  It  should  be  point¬ 
ed  out  that  the  Weber  number  defined  by  Leib  and  Goldstein 
is  the  inverse  of  that  used  in  this  work.  Lin  and  Lian*  showed 
that  absolute  instability  is  enhanced  by  the  presence  of  am¬ 
bient  fluid  in  a  sense  that  it  lowers  the  critical  Weber  num¬ 
ber.  They  proposed  that  absolute  instability  corresponds 
physically  to  the  dripping  mode  of  instability.  However, 
Tam  and  Hu9  did  not  find  absolute  instability  in  a  compress¬ 
ible  jet  emanating  at  high  supersonic  speeds  into  another 
compressible  fluid.  They  found  only  convective  instability. 
It  appears  that  the  compressibility  may  have  a  significant 
effect  on  the  onset  of  absolute  instability.  However,  while 
Tam  and  Hu  considered  only  the  case  of  a  high  supersonic  jet 
without  surface  tension,  the  previous  workers  neglected  the 
efTect  of  compressibility.  Thus  the  coupled  effects  of  surface 
tension  and  compressibility  remain  unclear.  These  effects 
are  very  important  in  many  practical  applications  involving 
a  high-speed  liquid  jet  emanating  into  a  gas  media.  For  ex¬ 
ample,  the  fuel  jet  in  certain  internal  combustion  engines 
may  reach  transonic  and  even  supersonic  speeds.  Here,  we 
investigate  the  effect  of  compressibility  on  the  instability  of  a 
liquid  jet  emanating  at  subsonic  as  well  as  supersonic  speeds 
into  an  ambient  fluid. 

II.  STABILITY  ANALYSIS 

Consider  the  onset  of  instability  of  a  jet  of  inviscid  com¬ 
pressible  fluid  emanating  from  a  circular  cylindrical  nozzle 


of  radius  a  into  an  unbounded  inviscid  compressible  fluid. 
The  conservation  equations  of  momentum,  mass,  energy, 
and  the  equation  of  state  are,  respectively, 

p.ld,  +Vf*)Vf«  -VP„  (1) 

0,  (2) 

id,  +  v,.v)e(  -  pm  +  y,'V)(prl)  =  o, 

and 

Pt=PM^)  (1=1,2), 

where  the  subscripts  /  =  1,2  denote,  respectively,  the  jet  fluid 
and  the  ambient  fluid,  p  is  density,  t  is  time,  V  is  the  velocity 
vector,  P is  the  dynamic  pressure,  e  is  the  internal  energy  per 
unit  mass,  and  T  is  temperature.  Note  that  the  flow  is  as¬ 
sumed  to  be  adiabatic.  The  basic  state  that  satisfies  Eqs.  ( 1 ) 
and  (2)  exactly  in  the  absence  of  gravity  is  given  by 

V,  =s  (UX*VX%WX)  =  (0,0, JF0), 

p]  =  const,  Px  =  const,  (0<R<fl), 

V,  =  (U2y^W2)  =  (0,0,0),  p2  =  const, 
P2=Pi—S/a  ( a<R<cc ), 

where  the  overbar  denotes  the  basic  state  quantities,  U,  V\ 
and  W  are,  respectively,  the  radial,  azimuthal,  and  axial 
components  of  velocity  in  the  cylindrical  coordinate  system 
(/?,#, Z),  W0  is  the  constant  jet  speed,  and  S  is  the  surface 
tension.  Z  is  positive  in  the  flow  direction. 

We  investigate  the  stability  of  this  basic  state  with  re¬ 
spect  to  the  disturbance  defined  by 

v,  =  v,  +  v;,  p,=P'+p;9  P,  =  P,  +  PU  r  =  a  +  d . 

(3) 

where  primes  are  used  to  designate  the  perturbations  of  the 
variables.  Substituting  Eq.  ( 3 )  into  Eqs.  ( 1 )  and  ( 2 ) ,  retain¬ 
ing  only  linear  terms,  we  have 

id,+6XtW0dz)\:=  -  (l/fr)V/\\  (4) 

(d,+6uW0dz)p:=  (5) 

where  <5,,  is  1  or  0  depending  on  if  /  =  1  or  iV  1  Taking  the 
divergence  of  (4)  and  then  using  Eq.  (5),  we  have 

id,  +6t,iv0dz)2p;  =  v2p;.  (6) 

Assuming  the  change  of  pressure  with  respect  to  density  is 

an  isentropic  process,  expanding  P,  [p,*s(prT, )  ]  about  p 
we  have  P\  =  (dP,/dp,  )$p '  =  c;p',  where  s  is  the  entropy. 
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Thus  Eq.  (6)  can  be  written  as 

(<?,  +  6U  W0  dz  )2P;  =  cfv2p;.  (7) 

With  respect  to  the  moving  fluids,  (7)  is  the  wave  equation 
with  the  speed  of  propagation  of  pressure  fluctuation  c, .  The 
speed  of  sound  c,  depends  on  the  specific  form  of  the  equa¬ 
tion  of  state  and  the  fluid  properties.  It  follows  from  the 
energy  equation  that 


dt  dt 


dh,  _  d(e,  +  P,/p, )  _  _  , dP, 
dt  dt  Pi  dt 

where  (d /dt)  stands  for  the  substantial  derivative  and 
h,  [  Tt>  st  (p,,  T{ )  ]  is  the  enthalpy.  Hence,  for  the  isentropic 
process  of  sound  propagation, 


f dhjdT ,\ 

),  p , 

V  de,  /dT, ) 

For  the  ideal  gas  surrounding  the  jet  P2  =  Rp2T2 ♦  dh2/dT2 
=  Cp,  de2/8T2  =  Cv9  and  thus  c2  =  ( yRT2)]/ 2,  where  R  is 
the  ideal  gas  constant,  and  y  =  Cp/Cv  is  the  ratio  of  the 
constant  pressure  specific  heat  Cp  to  the  constant  volume 
specific  heat  C, .  For  the  jet  liquid,  Px[(dhx/8Tx)/ 
(8ex/8Tx)]  =  K  is  the  adiabatic  bulk  modulus  of  elasticity, 
and  c,  =  (K /px)w2.  The  corresponding  linearized  bound¬ 
ary  conditions  for  Eq.  ( 7 )  at  the  interface  r  —  a  -f  d  are  the 
dynamic  force  balance 

P\  —  Pi  =  S(a~2  —  a  ~  2  dHe  —  dzz  )d  (8) 

and  the  kinematic  conditions 


(8t  +  6uWK%)2d**  —  p,  'dRP[,  (9) 

Note  that  P\  in  Eqs.  (8)  and  (9)  are  to  be  evaluated  at  r  =  a. 

The  normal  mode  solution  of  the  differential  system 
equations  (7)-(9)  will  be  sought  in  the  form 


P'2 

= 

p,H':oP:(r) 

af 

•exp[/(Az  +  nd  —  cor)}. 


(10) 


where  (r.z,/)  and  r  are  dimensionless  variables  defined  by 
(r9z,f)  =  (R.Z,d)a~i  and  T=t(W{/a)\ 

and  ( A,n)  and  co  are  the  dimensionless  wave-number  vector 
and  wave  frequency,  respectively. 

Substitution  of  Eq.  (10)  into  Eqs.  (7)-(9)  yields  the 
governing  equation  for  the  dimensionless  perturbation  pres¬ 
sure  amplitudes 


{D-  +  r~ 

'D-  [nV-J  +  k2 

-  M](k  —  a):]  }p, 

=  0. 

(11) 

{D-  +  r~ 

'D  —  |*7  “ 1  +  k  ‘ 

—  M  :w:  ]  }p2  =  0, 

(12) 

and  the  corresponding  boundary  conditions  at  r  = 

1, 

P:~P 

,+  »',(*  = +  n:- 

II 

O 

(13) 

o 

II 

1 

£ 

I 

(14) 

Qco:f- 

© 

II 

^  f*4 

(15) 

where  primes  denote  differentiation  with  r,  and 


K  = — Q  =  h  Mx  =  — ,  3/=-^. 
P\WoO  px  c,  c2 


Note  that  Wg  is  the  Weber  number,  Mx  is  the  Mach  number 
of  the  jet  and  M  is  the  Mach  number  of  the  ambient  fluid 
with  respect  to  a  frame  of  reference  moving  with  the  jet  at  a 
speed  WQ.  The  bounded  solution  of  Eqs.  (11)  and  (12)  are 


Pl=AInair),  Ax=[k2-M](k-co)2}'/2f 
p2  =-  BK„  (A2r)f  A 2=(k2-MW)w\ 


where  A  and  B  are  the  integration  constants  and  /„  and  Kn 
are,  respectively,  the  modified  Bessel  function  of  the  first  and 
second  kinds.  Substitution  of  these  solutions  into  Eqs.  (14) 
and  (15)  yields 

A={k-co)2f/rn(Ax)Ax, 

B=Qco2f/K'n(A2)A2 , 


where  primes  now  denote  the  differentiation  d  /d{Atr).  It 
follows  from  Eq.  (13)  that 


D  =  (k-co)2IjAx)/Air„(Al) 

-QorKM/WA*  2) 

-  We(k 2  +  nr  —  1)  =0.  (16) 


For  a  given  set  of  four  flow  parameters  Wt,  Q%  Mx,  and  M , 
this  is  the  secular  equation  the  solution  of  which  gives  the 
characteristic  wave  frequency  co  and  the  wave  number  k. 

For  the  special  case  of  an  incompressible  jet  in  a  vacu¬ 
um,  Q  =  0,  c,  —  oc ,  c2->  0,  thus  A ,  —  A,  A2  —  oc  and  Eq  (16) 
is  reduced  to,  if  co  remains  finite, 

/*  ( A)(A  —  co)2/l  ’n  (A)  +A(1  -  k2  -n2)Wt  =  0.  (17) 

This  is  the  characteristic  equation  used  by  Keller  et  a! 2  and 
Leib  and  Goldstein1  in  their  studies  of  the  spatial  instability 
of  the  Rayleigh  jet.  For  the  special  case  of  an  incompressible 
jet  emanating  into  an  incompressible  fluid, 
c,  —  oc,  Mx  —0,  thus  Ax  =  A2  =  A,  and  Eq.  ( 16)  is  re¬ 
duced  to,  if  co  remains  finite, 

(k  -  co)2In(k)/I  ’n(k)  -  Qco:Kn{k)/K  ’n(k) 

1  )WC  =0.  (18) 

This  equation  agrees  with  Eq.  ( 1 )  of  Lin  and  Lian*  w  hen  the 
Reynolds  number  of  their  equation  is  allowed  to  approach 
infinity  and  their  co  and  A  are,  respectively,  replaced  by  —  ico 
and  —  A.  This  change  of  variable  is  necessary  because  of  the 
different  positive  directions  and  the  different  definitions  of 
the  wave  frequency  that  are  used  in  these  two  papers  When 
the  surface  tension  is  put  to  zero,  Eq.  (18)  reduces  to  the 
characteristic  equation  corresponding  to  that  of  Tam  et  ai ,w 
except  their  equation  is  in  terms  of  the  Hankel  function  and 
the  Bessel  function  of  the  first  kind.  It  is  known  that  the 
axisymmetric  disturbance  grows  the  fastest  in  an  incom¬ 
pressible  jet  in  the  absence^  or  presence10  of  the  ambient 
fluid.  Here,  we  assume  that  the  same  is  true  for  a  compress¬ 
ible  jet.  Hence  we  treat  the  case  of  n  —  0  in  Eq  (16)  Here, 
both  co  ~  co,  4-  ico,  and  A  =  A,  -f  /A,  are  complex  Their 
imaginary  parts  co,  and  A,  give,  respectively,  the  temporal 
and  the  spatial  growth  rates  Their  real  parts  are  the  wa\e 
frequency  and  the  wave  number,  respectively. 
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III.  RESULTS 

The  numerical  accuracy,  as  well  as  the  correctness  of  the 
computer  program  that  solves  Eq.  ( 1 6 ) ,  were  ascertained  by 
comparing  some  known  results  of  published  works3*8 10  with 
our  results  for  the  corresponding  special  cases  These  known 
results  will  not  be  reproduced  here,  however.  Figure  1  delin¬ 
eates  the  solution  curves  of  (16)  for  the  case  of  an  incom¬ 
pressible  liquid  jet  emanating  into  a  gas  of  a  density  that  is 
1.3%  of  the  liquid  density.  This  density  ratio  corresponds  to 
a  water-air  system  under  1  atm.  The  Mach  number  of  the 
gas  relative  to  the  reference  frame  moving  with  the  liquid  jet 
is  supersonic.  Note  the  existence  of  the  two  saddle  points  C ' 
and  C,  respectively,  associating  with  the  two  sets  of  zero- 
temporal  growth-rate  curves  (A  'JJ')  and  (A>B).  The  con¬ 
stant  <o,  curves  initiate  from  the  upper  or  lower  k,  —  k, 
plane  with  finite  positive  values  of  to,.  The  values  of  to,  de¬ 
crease  toward  zero  before  they  reach  (A  \B ')  or  (AJi),  the 
curves  set  of  to,  =  0.  Here,  A  and  A  ’  will  be  called  the  down¬ 
stream  branches,  because  to,  increases  with  k,  near  these 
curves  and  thus  the  group  velocities  of  the  disturbances  ai  c 
in  the  downstream  direction.  Similarly,  the  disturbances 
corresponding  to  the  upstream  branches  B  and  B '  have  up¬ 
stream  propagating  group  velocities.  Observe  that  the  saddle 
point  C  has  the  value  of to,  lying  between  0.9603  and  0.9650. 
The  value  of  to,  at  C '  must  be  between  0  7443  and  0.75 1 6.  At 
these  saddle  points  D  =  dD  / dk ,  —  0  and  the  group  veloci  y 
is  zero.  Moreover,  co,  >  0  at  these  saddle  points.  The  appear¬ 
ances  of  the  saddle  pinch  point  signifies  that  the  Fourier 
integral  representation  of  the  disturbance  has  encountered  a 
singularity  along  its  path  of  integration.  Hence  the  integra¬ 
tion  path  cannot  be  analytically  deformed  long  before  co,  =0 
is  reached.  Consequently,  the  disturbance  will  grow  both 
temporj1’)  and  spatially  at  the  onset  of  absolute  instability 


FIG  1  Absolute  instabihtv  of  .ncompressible  jet  with  two  "addle  points 


as  explained  by  Bnggs’  and  Bers/  Note  that  t  lies  in  the 
region  of  k,  <  1  but  C  lies  in  the  region  of  k,  >  1.  The  absolute 
instabilities  associated  with  C '  and  C  will  be  termed  the  first 
and  ihe  second  kind,  respectively.  The  physical  significance 
of  this  will  be  explained  shortly.  When  M  is  reduced  to  sub¬ 
sonic  values  such  that  0  <  M  <  1 ,  with  the  rest  of  the  param¬ 
eters  specified  in  Fig.  1  fixed,  we  found  that  the  saddle  points 
C #  and  C  move  to  new  positions  but  remain  in  the  respective 
regions  of  k,  <  1  and  k,  >  1.  The  situation  remains  the  same 
even  when  A/—0.  However,  Lin  and  Lian8  found  that,  when 
M]  —  M  =  0,  there  is  only  one  saddle  point  belonging  to  the 
first  kind.  The  existence  of  the  second  saddle  point  in  this 
work  as  A/->0  therefore  consistutes  an  apparent  inconsis¬ 
tency  with  the  work  of  Lin  and  Lian.8  It  turns  out  that  Eq. 
(18)  used  by  Lin  and  Lian  and  Eq.  (17)  used  by  Leib  and 
Goldstein  are  both  singular  limits  of  Eq.  (16)  as  A/,  and 
A/-*  0.  The  solution  that  gives  rise  to  the  second  saddle  point 
cannot  be  recovered  from  Eqs.  (17)  and  (18),  because  Mxto 
and  Mto  remain  finite  as  A/,  and  A/— 0  in  Eq.  ( 16). 

Lin  and  Lian  also  showed  that  the  absolute  instability 
can  be  reduced  to  convective  instability  by  reducing  the  We¬ 
ber  number  if  both  the  jet  and  the  ambient  fluid  are  incom¬ 
pressible.  This  is  not  the  case  when  the  ambient  fluid  is  com¬ 
pressible  with  the  flow  parameters  specified  in  Fig.  1.  The 
values  of  flow  parameters  in  Fig.  2  remain  the  same  as  those 
in  Fig.  1,  except  that  We  is  reduced  from  1.5  to  0.1.  As  can  be 
seen  from  these  figures,  the  saddle  point  C'  can  indeed  be 
eliminated  by  reducing  W(  as  in  the  incompressible  case,  but 
the  new  saddle  point  C  remains.  Even  w  hen  We  =  0,  the  new 
saddle  point  still  remains.  With  A/=  1.74,  Wt  =0(10'') 
for  most  common  liquids,  and  thus  the  absolute  instability  of 
the  second  kind  should  be  easily  observable.  The  mecha¬ 
nisms  of  the  two  kinds  of  absolute  instabilities  are  different. 
A  discussion  on  these  two  different  mechanisms  will  be  pre¬ 
sented  near  the  end  of  this  section 


CO  10  20  3 C 
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FIG  2  Suppression  of  the  saddle  po  .it  of  th«.  first  kind  in  an  incompressible 
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An  extensive  computation  over  a  wide  range  of  flow 
parameters  shows  that  the  absolute  instability  of  the  second 
kind  can  be  suppressed  by  increasing  the  compressibility  of 
the  jet  while  the  rest  of  the  parameters  are  kept  constant. 
Figures  3  and  4  provide  such  an  example.  For  the  flow  pa¬ 
rameters  specified  in  Fig.  3,  there  are  both  kinds  of  saddle 
points.  The  Saddle  point  of  the  second  kind  is  eliminated  in 
Fig.  4  by  further  increasing  Af,  to  0.52  from  0.51  in  Fig.  3. 
However,  fairly  extensive  computations  show  that  the  elimi¬ 
nation  of  the  saddle  point  of  the  second  kind  can  be  achieved 
only  when  M  >  1.  Note  that  the  saddle  point  of  the  first  kind 
still  remains.  The  saddle  point  of  the  first  kind  persists  even  if 
one  continues  to  increase  M]  to  a  high  supersonic  value.  This 
is  demonstrated  in  Fig.  5.  It  should  be  pointed  out  that  the 
amplification  curve  for  the  convectively  unstable  second 
mode  appearing  in  Fig.  4  still  remains  when  the  flow  param¬ 
eters  are  changed  to  that  of  Fig.  5.  It  does  not  appear  in  Fig. 
5,  because  it  lies  outside  of  the  right  edge  of  this  figure. 

When  the  surface  tension  of  the  jet  characterized  in  Fig. 
5  is  reduced  such  that  Wt  is  decreased  from  0.8  to  0.6,  the 
absolute  instability  is  completely  eliminated,  and  the  jet  be¬ 
comes  convectively  unstable  as  shown  in  Fig.  6.  The  a>t  =  0 
curve  in  Fig.  6  can  be  approached  along  constant  cor  curves 
originated  with  ot  >  0  from  the  upper  k,-kt  plane.4  ?  Hence 
the  curve  cd,  —  0  in  Fig.  6  gives  the  spatial  amplification  rate 
of  the  convectively  unstable  disturbances  Further  reduction 
of  does  not  reintroduce  the  absolute  instability.  This  is 
consistent  with  the  results  of  Tam  and  Hu  who  found  no 
absolute  instability  at  high  supersonic  speeds  of  the  jet  with¬ 
out  surface  tension. 

In  order  to  understand  the  physical  mechanisms  of  the 
absolute  instability  and  of  convective  instability,  we  obtain 
the  eigenvectors  (p\,p2f)  corresponding  to  the  eigenvalues 
near  the  two  types  of  saddle  points.  The  pressure  exerted  by 
the  jet  on  its  interface  corresponding  to  the  absolute  mstabil- 


FIG  3  Absolute  instability  of  a  supersonic  jet 


FIG.  4.  Suppression  of  the  saddle  point  of  the  second  kind 


ity  of  the  first  kind  depicted  in  Fig.  4  is  plotted  in  Fig.  7, 
together  with  p2  and  the  spatially  growing  interfacial  dis¬ 
placement  /  The  pressure  distribution  is  nearly  180°  out  of 
phase  with  the  displacement.  Thus  ( p}  —  p2)  stabilizes  the 
jet  by  acting  against  necking  and  bulging,  respectively,  at  the 
wave  troughs  and  wave  peaks.  On  the  contrary ,  the  surface 
tension  force  tends  to  facilitate  capillary  pinching  along  the 
radial  direction.  This  can  be  seen  from  Eq.  ( 16).  Note  that 
the  —  k  2  We  and  the  Wc  terms  in  Eq.  (16)  have  opposite 
signs.  It  is  seen  from  Eq.  (8)  that  the  —  k  2  We  term  arises 
from  the  interfacial  curvature  along  the  axial  direction,  and 
the  W'  term  arises  from  the  curvature  along  the  circumfer¬ 
ential  direction.  While  the  surface  tension  uses  the  curvature 


FIG  5  Absolute  instability  of  supersonic  jet 
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FIG  6  Convective  instability  of  supersonic  jet 


along  the  axial  direction  to  restore  the  interface  to  its  unper¬ 
turbed  position,  it  exploits  the  curvature  along  the  circum¬ 
ferential  direction  to  pinch  the  jet.  For  the  absolute  instabil¬ 
ity  of  the  first  kind  k  <  1  near  the  saddle  point,  and  thus  the 
—  k  2  We  term  is  dominated  by  the  W9  term.  Hence  the  jet 
becomes  unstable  due  to  the  capillary  pinching.  For  the  ab¬ 
solute  instability  associated  with  the  saddle  point  of  the  sec¬ 
ond  kind,  k  >  1  and  thus  surface  tension  is  a  stabilizing  force. 
The  instability  must  be  caused  by  the  unfavorable  pressure 
distribution.  The  pressure  distribution  along  the  interface 
corresponding  to  the  absolute  instability  of  the  second  kind 
depicted  in  Fig.  2  is  presented  in  Fig.  8.  Note  that  (/?,  —  p2) 
in  Fig.  8  tends  to  steepen  the  slope  of  the  wave  front.  Appar¬ 


F1G  7  Absolute  instabilit)  of  the  first  kind 


FIG  8.  Absolute  instability  of  the  second  kind 


ently,  the  stabilizing  surface  tension  is  not  sufficiently  large 
to  overcome  the  pressure  force. 

It  is  seen  from  Fig.  7  that  the  pressure  in  the  jet  opposes 
the  first  kind  of  absolute  instability.  The  mean  pressure  in 
the  jet  is  inversely  proportional  to  Q  It  follow  s  that  a  larger 
value  of  Q  may  promote  the  absolute  instability  of  the  first 
kind.  Figures  6  and  9  demonstrate  this  point  clearly.  The 
convective  instability  depicted  in  Fig.  6  changes  over  to  the 
absolute  instability  portrayed  in  Fig.  9  simply  because  of  the 
increasein  Q  from  0.0013  to  0.08.  It  was  also  found  that  the 
absolute  instability  of  the  second  kind  has  the  same  depend¬ 
ence  on  Q. 


FIG  9  Absolute  instability  of  the  first  kind 
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IV.  CONCLUSION 

In  the  four-dimensional  hyperspace  spanned  by  Qy 
A/,,  and  My  there  exists  a  surface  that  separates  the  space  into 
regions  of  convective  and  absolute  instabilities.  A  simple  de¬ 
lineation  of  this  surface  cannot  be  made  without  a  complete 
mapping  of  this  parameter  space.  The  stability  characteris¬ 
tics  of  a  compressible  jet  encountered  in  any  specific  applica¬ 
tions  have  to  be  determined  from  the  solutions  of  Eq.  (16) 
for  the  specific  flow  parameters  encountered.  However, 
based  on  the  fairly  extensive  computation  carried  out  by  us, 
we  may  make  the  following  points.  Two  kinds  of  absolute 
instability  were  found.  The  first  kind  is  caused  by  the  capil¬ 
lary  pinching,  but  the  second  is  due  to  the  interfacial  pres¬ 
sure  fluctuation  alone.  A  reduction  of  surface  tension  may 
eliminate  the  first  kind,  but  cannot  remove  the  second  kind. 
The  compressibility  of  the  jet  may  be  used  to  deviate  the 
absolute  instability  of  the  second  kind  if  M>  1.  The  com¬ 
pressibility  of  the  ambient  fluid  promotes  the  absolute  insta¬ 
bility.  The  absolute  instability  of  both  kinds  can  be  removed 
simultaneously  by  use  of  lower  surface  tension  and  larger  jet 
compressibility.  Whether  the  instability  of  a  jet  is  absolute  or 
convective  depends  on  the  dispersion  property  dictated  by 
the  particular  parameter  range  encountered  by  the  jet.  When 
the  spatially  growing  disturbance  is  capable  of  escaping  with 
a  finite  group  velocity  from  a  point  in  space  the  jet  is  convec- 
tively  unstable,  otherwise  it  is  absolutely  unstable.  As  far  as 
the  physical  mechanism  is  concerned,  both  convective  and 
absolute  instabilities  are  either  due  to  capillary  pinching  or 
the  pressure  fluctuation. 
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Absolute  and  convective  instability  of  a  viscous  liquid  jet  surrounded 
by  a  viscous  gas  in  a  vertical  pipe  C  APpGNP<x  V  ) 
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The  absolute  and  convective  instability  of  a  viscous  liquid  jet  emanating  into  a  viscous  gas  in 
a  vertical  pipe  is  analyzed  in  a  parameter  space  spanned  by  the  Reynolds  number,  the 
Froude  number,  the  Weber  number,  the  viscosity  ratio,  the  density  r ?*'*  and  the  diameter 
ratio.  The  numerical  results  of  the  analysis  are  used  to  demonstrate  that  reduction  in 
gravity  tends  to  enhance  the  Rayleigh  mode  of  convective  instability  which  leads  to  the 
breakup  of  a  liquid  jet  into  drops  of  diameters  comparable  with  the  jet  diameter.  On 
the  contrary,  the  Taylor  mode  of  convective  instability  that  leads  to  atomization  is  retarded 
at  reduced  gravity.  The  Rayleigh  mode  becomes  absolutely  unstable  when  the  Reynolds 
number  exceeds  a  critical  value  for  a  given  set  of  the  rest  of  the  relevant  parameters.  The 
domain  of  absolute  instability  is  significantly  enlarged  when  the  effect  of  gas  viscosity 
is  not  neglected. 


The  stability  of  core-annular  flows  has  been  investi¬ 
gated  extensively  in  the  context  of  lubricated  pipelining1  as 
well  as  in  the  context  of  liquid  jet  atomization.2  Tempo¬ 
rally  growing  disturbances  were  considered  in  the  pub¬ 
lished  works.  In  the  present  work,  we  investigate  the  sta¬ 
bility  of  core-annular  flow  with  respect  to  spatially  growing 
disturbances.  Spatially  growing  disturbances  are  especially 
relevant  to  the  atomization  process  that  utilizes  the  ampli¬ 
fication  of  disturbances  as  they  are  convected  downstream 
to  generate  small  droplets.  The  atomization  process  must 
be  carried  out  outside  of  the  Rayleigh  regime3  of  jet 
breakup  by  capillary  pinching.  It  must  also  avoid  absolute 
instability.4,5  The  critical  Weber  number  above  which  the 
flow  is  absolutely  unstable  is  obtained  as  a  function  of  the 
Reynolds  number  for  given  sets  of  the  rest  of  relevant  flow 
parameters.  Outside  of  the  region  of  absolute  instability  in 
the  parameter  space,  it  is  shown  that  there  are  two  physi¬ 
cally  distinct  modes  of  convective  instability. 

Consider  the  stability  of  a  cylindrical  liquid  jet  of  den¬ 
sity  p]  flowing  vertically  downward  in  the  direction  of 
gravitational  acceleration  g.  The  jet  is  surrounded  by  a  gas 
of  density  p2  occupying  the  annular  space  between  the  jet 
and  the  inner  wall  of  a  vertical  pipe  of  radius  R2  The  liquid 
and  gas  flows  have  the  same  pressure  gradient  in  the  jet 
flow  direction.  Hence  the  pressure  drop  across  the  liquid- 
gas  interface  at  any  cross  section  is  constant  along  the  axial 
direction.  This  pressure  difference  is  balanced  exactly  by 
the  surface  tension  force  per  unit  surface  area  of  the  liquid 
cylinder  of  constant  radius  R]f  which  is  also  the  radius  of 
curvature.  The  velocity  distribution  of  the  basic  flow  that 
satisfies  exactly  the  Navier-Stokes  equations  is  given  by 
Eq.  ( 1 )  of  Lin  and  Ibrahim.2  The  linear  stability  analysis 
of  this  basic  flow  with  respect  to  the  normal  mode  distur¬ 
bance  yields  the  characteristic  equation  given  by  Eq.  (21 ) 
of  Lin  and  Ibrahim.2  The  normal  mode  disturbances  in¬ 
volve  an  exponential  factor  exp (//cz +&>/),  where  t  is  time 
and  z  is  the  axial  distance  measured  in  the  opposite  direc¬ 


tion  of  the  liquid  flow.  The  wave  number  k  was  treated  as 
real  in  the  previous  analysis2  of  temporally  growing  distur¬ 
bances.  In  this  work,  both  k  and  o  are  complex,  i.e., 
k  =  kr+ikn  and  iYu,.  Hence  ^>0  gives  the  spatial 

growth  rate  of  disturbances  convected  as  a  group  in  the 
downstream  direction.  Here,  <ur>0  gives  the  temporal 
growth  rate  of  unstable  disturbances,  kr  is  the  real  wave 
number  that  is  inversely  proportional  to  the  wavelength, 
and  o)i  is  the  natural  frequency  of  oscillation  of  the  normal 
mode  disturbance.  The  group  velocity  of  the  disturbance  is 
given  by  —do)fdkr 

To  determine  the  spatial  amplification  rate  of  convec- 
tively  unstable  disturbances  for  a  given  set  of  parameters, 
we  put  cor=  0  and  assign  different  values  to  co(  in  the  char¬ 
acteristic  equation.  Then,  kr  and  k,  are  solved  from  the  real 
and  imaginary  parts  of  the  characteristic  equation  by  use  of 
a  modified  Powell’s  hybrid  method  with  the  subroutines 
CGEDI  and  DNSQ  in  the  SLATEC  library*.  The  set  of  points 
(kn  kr)  for  various  values  of  o,  determine  the  spatial  am¬ 
plification  curve  <ur=0.  The  number  of  terms  retained  in 
the  series  solutions  are  varied  and  the  resulting  character¬ 
istic  equations  are  solved  for  the  eigenvalues.  The  charac¬ 
teristic  values  corresponding  to  different  orders  of  the  char¬ 
acteristic  determinant  are  compared  to  ascertain  the 
desired  numerical  accuracy.  It  is  found  that  it  is  sufficient 
to  retain  only  three  and  five  terms,  respectively,  in  the 
liquid  and  gas  phases  for  four  decimal  point  accuracy  for 
the  given  parameters. 

The  numerical  accuracy,  as  well  as  the  correctness  of 
the  computer  program  that  solves  the  characteristic  equa¬ 
tion,  were  ascertained  by  comparing  some  known  results  of 
published  works4  5  with  our  results  for  the  corresponding 
special  cases.  These  known  results  will  not  be  reproduced 
here,  however.  The  relevant  parameters  are  the  Weber, 
Reynolds,  and  Froude  numbers,  the  density,  viscosity  and 
radius  ratios  defined,  respectively,  by  We=S/p, 

R  c  =  p\W'0R]/p],  W\fgRvp2/pv  p2/pv  R2/Rv  where  W0 
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FIG.  1.  The  existence  of  the  absolute  instability. 

is  the  maximum  velocity  in  the  liquid  jet,  fi]  and  \i2  are> 
respectively,  the  dynamic  viscosities  of  the  jet  liquid  and 
the  surrounding  gas. 

The  results  on  absolute  instability  that  were  not  found 
in  the  temporal  analyses  will  be  reported  first.  For  a  suffi¬ 
ciently  high  Weber  number,  given  the  rest  of  relevant  pa¬ 
rameters,  we  obtain  two  branches  of  spatial  amplification 
curves  cjr=0.  An  example  is  shown  in  Fig.  1,  where  the 
two  branches  are  designated  with  A  and  B.  The  values  of 
Oj  increase  along  A  and  decrease  along  B.  Thus  the  group 
velocities  of  disturbances  associated  with  branches  A  and  B 
propagate,  respectively,  in  the  downstream  and  upstream 
directions.  Notice  the  appearance  of  a  saddle  point  in  this 
figure.  The  constant  curves  initiated  from  the  upper  or 
lower  kr-k,  plane  with  finite  positive  values  of  <or  delineated 
the  manner  with  which  the  saddle  point  is  formed.  Observe 
that  the  saddle  point  has  the  value  of  wl  lying  between  0.73 
and  0.74.  At  the  saddle  point  D=dD/dkr=0 ,  where  D 
stands  for  the  characteristic  determinant,  i.e.,  the  left  side 
of  the  characteristic  equation.  At  the  saddle  point,  the 
group  velocity  is  zero,  and  G>r>0.  The  appearance  of  the 
saddle  point  signifies  that  the  Fourier  integral  representa¬ 
tion  of  the  disturbance  has  encountered  a  pinch  point  sin¬ 
gularity  along  its  path  of  integration.  Hence  the  integration 
path  cannot  be  analytically  deformed  long  before  <ur=0  is 
reached.  Consequently,  the  disturbances  will  grow  both 
temporally  and  spatially  at  the  onset  of  absolute  instability 
as  explained  by  Briggs  and  Bers.6  When  the  value  of  We  is 
gradually  decreased  with  the  rest  of  the  parameters  in  Fig. 

1  kept  constant,  the  saddle  point  moves  gradually  toward 
the  downstream  branch  and  eventually  disappears.  Thus 


FIG.  2.  Critical  Weber  number  of  absolute  instability 

there  exists  a  critical  Weber  number  above  which  the  flow 
is  absolutely  unstable,  and  below  which  it  is  convectively 
unstable  for  a  given  set  of  parameters.  The  critical  Weber 
number  is  found  to  be  only  weakly  dependent  on  Fr,  A\ 
and  /.  The  critical  Reynolds  number  as  a  function  of  We 
for  2=0.0013,  A=0.018,  Fr=  1000,  and  /=  10  is  depicted 
as  curve  B  in  Fig.  2.  Curve  A  in  this  figure  is  the  corre¬ 
sponding  curve  for  the  case  of  inviscid  surrounding  gas  of 
infinite  extent,5,  i.e.,  for  the  case  of  jV=0,  /=  oc,  and 
2=0.0013.  It  is  observed  that  the  domain  of  absolute  in¬ 
stability  is  significantly  enlarged. 

Figure  3  delineates  the  dependence  of  the  spatial  am¬ 
plification  curves,  o>r=0,  on  the  Weber  number  for  the 
given  set  of  parameters  specified  in  the  figure.  This  set  of 
parameters  falls  inside  the  region  of  convective  instability. 
All  of  the  amplification  curves  can  be  reached  along  con¬ 
stant  a>i  curves  initiating  from  the  lower  krk,  plane  with 
finite  positive  values  of  wr  Moreover,  the  group  velocity  of 
the  disturbance,  - ~da)/dkn  is  negative  along  each  of  the 
amplification  curves.  Hence  these  curves  give  the  spatial 
amplification  rate  of  the  disturbances  as  they  are  swept 
downstream  as  a  group.  It  is  seen  that  the  amplification 
rate  increases  with  We,  indicating  that  surface  tension  is  a 
source  of  the  convective  instability.  It  is  also  seen  that  the 
wavelength  corresponding  to  the  most  amplified  distur¬ 
bance  is  of  the  same  order  of  magnitude  as  the  jet  radius,  as 
is  encountered  in  the  Rayleigh  jet.  These  resists  are  qual¬ 
itatively  the  same  as  those  found  in  the  temporal  theory  .2 
Note  that  the  values  of  N  and  Q  used  in  Fig.  3  correspond 
to  a  water  jet  in  air  at  room  temperature  under  one  atmo¬ 
spheric  pressure.  The  effects  of  the  rest  of  the  parameters 
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FIG.  3.  The  Rayleigh  mode  convective  instability. 

on  the  convectively  unstable  disturbances  remain  qualita¬ 
tively  the  same  as  those  for  the  case  of  temporally  growing 
disturbances.  Therefore  they  will  not  be  explained  here. 
However,  the  effect  of  the  Weber  number  of  Taylor  mode 
instability  at  microgravity  will  be  mentioned  here,  since  it 
was  not  pointed  out  previously.  When  the  parameters  in 
Fig.  3  are  kept  the  same  except  the  Froude  number  is 
reduced  to  2000,  the  jet  becomes  unstable  with  respect  to 
the  Taylor  mode.  Contrary  to  the  situation  in  Fig.  3,  the 
amplification  rate  decreases  with  We,  indicating  that  sur¬ 
face  tension  force  is  not  the  source  of  instability.  The  other 
noticeable  feature  is  that  the  wavelength  corresponding  to 
the  most  amplified  disturbance  for  this  set  of  parameters  is 
one  order  of  magnitude  smaller  than  the  jet  diameter,  as  is 
frequently  observed  in  atomization  processes.  The  physical 
mechanism  of  atomization  was  shown3  to  be  the  interfacial 
pressure  fluctuation  for  the  case  of  a  viscous  jet  in  an  in- 
viscid  gas.  Here,  with  the  gas  viscosity  included,  the  inter- 
facial  pressure  fluctuation  is  supplemented  by  the  viscous 
stress  fluctuation  to  enhance  the  instability.  This  enhance¬ 
ment  can  be  expected  to  be  significant  when  the  change  of 


velocity  gradient  across  the  interface  becomes  large.  The 
change  of  the  velocity  gradient  increases  with  the  decrease 
in  the  value  of  Fr  as  can  be  seen  from  the  expressions  of  the 
basic  state  velocity  distribution.  This  is  a  physical  explana¬ 
tion  of  why  the  Rayleigh  mode3  of  instability  depicted  in 
Fig.  3  changes  to  Taylor’s3  mode  of  atomization  instability 
at  reduced  gravity. 

The  liquid  core-gas  annular  flow  of  two  viscous  fluids 
is  shown  to  be  capable  of  exhibiting  either  absolute  or  con¬ 
vective  instability  depending  on  the  specific  values  of  flow’ 
parameters  involved.  The  critical  Weber  number  below 
which  the  flow  is  convectively  unstable  and  above  which 
the  flow  is  absolutely  unstable,  is  obtained  as  a  function  of 
the  Reynolds  number  for  the  given  remaining  parameters. 
The  gravity  is  found  to  have  little  effect  on  the  absolute/ 
convective  nature  of  instability.  There  are  two  modes  of 
convective  instability:  the  Rayleigh  mode  of  convective  in¬ 
stability  caused  by  capillary  pinching  and  the  Taylor  mode 
of  atomization  caused  by  interfacial  stress  fluctuation.  The 
Taylor  mode  instability  is  shown  to  change  over  to  the 
Rayleigh  mode  when  the  Froude  number  is  increased  be¬ 
yond  a  certain  value  with  the  rest  of  flow  parameters  fixed. 
The  practical  implication  of  this  is  that  the  atomization 
processes  widely  used  on  earth  in  a  given  parameter  range 
may  fail  to  produce  droplets  smaller  than  the  jet  thickness 
at  microgravity.  In  order  to  recover  the  atomization  pro¬ 
cess  at  microgravity,  one  may  reduce  the  width  of  the  an¬ 
nular  gas  layer  to  increase  the  interfacial  stress. 
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Abstract 


Current  fuel  injection  design  relies  upon  empirical  studies  to  optimize  injector 
technology.  The  fundamental  understanding  *he  relationship  between  the  fuel  pres¬ 
sure  input  to  an  injector  nozzle  and  the  resulting  spray  has  remained  vague.  An  experi¬ 
mental  study  was  undertaken  to  measure  the  spray  characteristics  of  a  diesel  spray  re¬ 
sulting  from  various  dynamic  pressure  inputs  to  a  nozzle. 

An  intermittent  liquid  spray  was  created  by  applying  successive  cycles  of  a  pres¬ 
sure  waveform  to  a  plain  hole  nozzle.  Drop  sizes  and  their  temporal  relationship  to  the 
pressure  waveform  were  recorded  at  various  locations  downstream  from  the  nozzle. 
The  results  are  presented  for  several  nozzle  pressure  waveforms. 
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Chapter  1 


Introduction 

In  order  to  understand  fuel  spray  formation  in  fuel  injection  systems,  this  work  has 
investigated  the  basic  spray  formation  of  a  simple,  intermittent  fuel  spray.  Several  dif¬ 
ferent  nozzle  pressure  waveforms  were  applied  to  a  single,  cylindrical  orifice.  The  re¬ 
sulting  spray  was  then  analyzed  using  Aerometric’s  Phase/Doppler  Particle  Analyzer, 

The  rise  and  fall  of  the  driving  pressure  in  the  nozzle  subjects  the  fluid  to  differ¬ 
ent  rates  of  acceleration  as  it  leaves  the  nozzle.  Thus,  if  different  pressure  pulses  are 
applied  to  the  nozzle,  there  should  be  some  consistent  pattern  of  drop  size  and  veloc¬ 
ity  distributions  since  the  spray  formation  is  a  dynamic  process  which  is  sensitive  to  the 
dynamic  pressure  input.  Comparison  of  different  pressure  waveforms  in  the  nozzle  and 
their  respective  spray  distributions  can  then  be  used  to  optimize  intermittent  spray 
technology. 

There  are  several  studies  of  intermittent  sprays  which  measure  the  time-resolved 
drop  size  and  velocity  distributions  [1-5].  These  studies  have  all  relied  upon  conven¬ 
tional  fuel  injection  pumps  and  injectors  that  are  difficult  to  use  for  a  systematic  evalu¬ 
ation  of  a  dynamic  pressure  input’s  influence  on  spray  evolution.  These  systems  were 
typically  designed  for  specific  commercial  applications. 
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For  example,  Sangeorzan  et  al  [3]  give  time-resolved  spray  distributions  produced 
by  a  conventional  diesel  injection  system.  The  re  is  no  systematic  investigation  of  the 
influence  of  a  dynamic  pressure  input  to  the  injector,  but  they  point  out  that  there  are 
cycle  to  cycle  variations  in  the  spray  produced.  They  indicate  that  these  cycle  to  cycle 
variations  are  probably  due  to  line  pressure  variations.  However,  the  basic  influence 
of  a  dynamic  pressure  input  on  spray  formation  has  still  been  left  vague.  The  line  pres¬ 
sure  variations  were  a  defect  in  their 
system,  not  the  object  of  their  study. 

Nonconventional  injection  sys¬ 
tems  have  also  been  used  for  injection 
studies.  For  example,  Sinnamon  et  al 
[11]  used  a  solenoid-operated  injector 
which  creates  nearly  square  nozzle 
pressure  waveforms  of  various  dura¬ 
tions  during  their  study  of  spray  trajec¬ 
tories  in  an  engine.  Ruiz  and  Chigier  [12]  used  the  same  type  of  solenoid-operated  in¬ 
jector  to  study  the  spray  angle.  Although  they  do  not  measure  the  influence  of  a  dy¬ 
namic  pressure  input  on  the  drop  size  and  velocity  distributions  directly,  they  do  study 
its  indirect  influence  through  observed  spray  angles  and  trajectories. 

Althought  the  influence  of  a  dynamic  pressure  input  has  not  been  studied  directly, 
its  importance  has  been  suggested.  It  is  clear  from  Gliken  [6]  and  Dolenc  [7]  that  the 
dynamic  pressure  input  to  the  nozzle  is  very  important  to  engine  efficiency.  In  their  dis¬ 
cussion  of  optimizing  fuel  injection  systems,  they  point  out  how  injection  flow  rate  is 
responsible  for  such  qualities  of  engine  performance  as  specific  fuel  consumption, 
noise,  nitrous  oxide  emissions,  and  smoke  emissions.  General  design  criteria  show  the 
trade-offs  of  having  different  injection  rates  at  different  times  in  the  injection  process 


Figure  1.1  Engine  performance  trade-offs  in  tailoring 
fuel  injection  rate  within  the  time  of  a  single  pulse  [7], 


2 


(sec  figure  1.1).  Since  the  pressure  differentia]  across  the  nozzle  orifice  dictates  the  in¬ 
jection  rate,  we  can  conclude  the  importance  of  a  dynamic  pressure  input  to  engine  ef¬ 
ficiency. 

There  has  also  been  the  practical  application  of  injection  rate  tailoring  to  conven¬ 
tional  fuel  injection  pump,  pipe,  and  nozzle  systems.  Computer  programs  [6,8]  have 
been  written  which  can  predict  time-resolved  injection  rates  for  specific  systems.  En¬ 
gineers  use  these  programs  to  optimize  hardware  components  according  to  the  perfor¬ 
mance  trade-offs  like  those  in  figure  1.1.  Two-stage  injectors  [9,10]  have  also  been  de¬ 
veloped  to  provide  low  injection  rates  at  the  start  of  injection  and  higi.  injection  rates 
during  the  remainder  of  the  injection.  The  result  of  two-stage  injection  is  lower  noise 
and  emissions  for  engines  using  them. 

However,  by  focusing  on  injection  rate  alone  and  not  spray  formation  as  a  result 
of  a  particular  injection  rate  schedule,  a  fundamental  basis  for  engine  efficiency  is  not 
established.  The  fundamental  change  taking  place  with  varying  injection  rates  is  a 
change  in  spray  characteristics  that  effect  the  combustion  process.  Further  investiga¬ 
tion  should  be  directed  toward  the  connection  between  injection  rate  and  spray  char¬ 
acteristics,  as  well  as  the  impact  of  different  sprays  on  combustion.  This  is  the  ratio¬ 
nale  for  the  present  study. 

Further  complicating  the  issue  of  atomization  in  fuel  injection  is  the  problem 
posed  by  the  small  injection  chambers  found  in  reciprocating  internal  combustion  en¬ 
gines.  Studies  of  spray  angles  and  trajectories  are  important  because  they  are  con¬ 
cerned  with  optimizing  the  spray  within  the  confines  of  small  chambers.  Likewise,  it 
is  important  to  study  the  dense  core  of  intermittent  sprays  which  tends  to  splash  off  the 
chamber  walls.  Previous  experimental  studies  [13,14]  of  the  core  have  tried  to  measure 
the  size  of  the  core.  Koo  and  Martin  [5]  have  shown  that  Aerometric's  PDPA  system 
can  measure  drop  size  distributions  in  dense  regions  near  the  core  with  some  success. 
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Characterization  of  the  core  in  relation  to  different  pressure  inputs  is  the  key  to  mini¬ 
mizing  the  core  size. 

Optimization  of  the  injection  spray  in  general,  has  been  done  empirically. 
Successful!  injection  systems  are  adjusted  to  give  the  best  practical  performance  with¬ 
out  regard  to  the  atomization  process.  The  next  step  is  to  consider  the  physical  process 
of  atomization  in  the  optimization  of  injection.  Numerical  models  will  someday  be  able 
to  predict  optimal  design  criteria. 

The  influence  of  a  dynamic  pressure  input  upon  spray  formation  is  being  analyti¬ 
cally  investigated  by  Lin,  Woods,  and  Zhou.  They  are  developing  a  numerical  model 
from  first  principles  to  compute  the  motion  of  a  liquid  through  a  plain  orifice,  into  the 
core  region,  and  on  to  the  outer  atomized  regions.  Validation  of  their  work  through 
experimental  results  like  those  presented  here  will  tie  together  all  of  the  factors  in¬ 
volved  in  atomizing  a  liquid  jet  through  a  cylindrical  orifice. 
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Chapter  2 


Experimental  Setup  and  Procedure 


2.1  Experimental  Setup 

The  experimental  apparatus  was  designed  to  produce  regular  pulses  of  high  pres¬ 
sure  liquid  at  the  nozzle;  record  the  time-resolved  pressure  waveforms  of  the  high  pres¬ 
sure  pulses;  capture  drop  size  and  velocity  information  about  the  resulting  spray;  and 
contain  the  spray  in  a  safe,  non-obtrusive  manner.  It  was  of  particular  importance  that 
the  apparatus  be  able  to  generate  regular  pulses  with  a  user  selectable  pressure  wave¬ 
form  in  order  to  automate  the  process. 

Conventional  fuel  injection  equipment  was  not  used  mainly  for  its  inflexibility  in 
producing  user  selectable  nozzle  pressure  waveforms.  In  addition,  the  nozzle  pressure 
waveforms  produced  by  previous  studies  [1-5]  of  time-resolved  spray  analysis  were  not 
smooth.  Typically,  there  are  high  frequency  pressure  fluctuations  superimposed  upon 
a  square  pressure  waveform.  Although  these  fluctuations  are  of  practical  importance, 
it  is  more  difficult  to  investigate  them  analytically.  Smooth,  low  frequency  nozzle  pres¬ 
sure  waveforms  are  more  suitable  for  the  longer  time  steps  in  numerical  models  that 
may  use  this  data. 

Thus,  a  custom  system  was  built  around  an  electrohydraulic  valve  that  controlled 
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fuel  flow  proportionally  to  an  electric 
signal  in  order  to  meet  the  need  for  a 
smooth  pressure  waveform.  The  system 
consistently  produced  smooth,  low  fre¬ 
quency  nozzle  pressure  waveforms  as 
desired. 

The  basic  system  involves  a  high 
pressure  pump  that  supplies  diesel  at 
nearly  constant  pressure  to  the  electro- 
hydraulic  valve.  A  signal  generator  modulates  the  flow  through  the  valve  to  the  nozzle. 
The  nozzle  pressure  is  measured  by  a  transducer  located  in  a  port  behind  the  cylindri¬ 
cal  nozzle  orifice  (see  figure  2.2).  A  two-axis  traverse  moves  the  nozzle  relative  to  the 
fixed  probe  volume,  changing  the  spatial  sampling  location  in  the  spray.  The  follow¬ 
ing  paragraphs  describe  the  components  in  further  detail. 

A  high  pressure,  air-driven  liquid  pump  (Haskel  Model  AW-35)  was  used  to  pro¬ 
vide  a  source  of  high  pressure  fuel  to  the  electrohydraulic  valve.  This  pump  utilizes  a 
reciprocating  piston  to  move  and  compress  the  fluid.  The  reciprocating  action  produces 
an  unwanted  fluctuation  in  the  supply  pressure.  A  hydraulic  accumulator  just  down¬ 
stream  of  the  pump  is  an  effective  approach  to  reducing  these  unwanted  fluctuations 
which  are  common  to  all  piston  pumps.  In  this  case,  the  expense  of  an  accumulator  was 
avoided,  since  the  pump  provided  a  long  period  of  constant  pressure  compared  to  the 
period  of  unsteadiness.  Those  injections  associated  with  the  supply  pressure  unsteadi¬ 
ness  could  be  rejected  by  later  processing  because  their  maximum  pressure  amplitude 
is  outside  a  selected  range.  That  is,  injections  were  generated  without  regard  to  the 
pump  pressure  fluctuations.  The  data  associated  with  the  periods  of  supply  pressure 
unsteadiness  are  characterized  by  pressure  pulses  of  low  maximum  pressure  amplitude. 


Figure  2.2  Details  of  the  nozzle  design. 


Processing  of  the  data  involved  selecting  data  for  which  the  maximum  pressure  ampli¬ 
tude  falls  within  the  936psi  to  1036psi  range.  Only  those  pulses  selected  were  used  for 
the  statistical  studies. 

Originally,  the  reciprocating  action  of  this  pump  was  meant  to  be  the  direct  source 
of  intermittent  pulses;  however,  the  pump’s  low  frequency  and  the  limited  control 
achievable  over  the  pressure  waveform  proved  inadequate.  Instead,  the  intermittent 
pulse  was  generated  by  an  electrical  signal  to  a  proportional  hydraulic  valve.  This  valve 
would  open  and  close  proportionally  to  an  electrical  signal  within  its  limited  frequency 
response  and  modulate  the  flow  of  liquid  to  the  nozzle. 

A  signal  generator  (Hewlett-Packard  Model  3312A)  was  used  directly  to  control 
the  electrohydraulic  valve.  Both  the  frequency  and  the  waveform  of  the  electrical  sig¬ 
nal  could  be  controlled  with  this  generator,  thereby  permitting  control  of  the  intermit¬ 
tent  spray  pulsation  frequency  and  nozzle  pressure  waveform. 

The  nozzle  pressure  was  recorded  with  a  digital  data  acquisition  system  (Keithley 
Model  S-500)  using  a  piezoelectric  pressure  transducer  (PCB  Piezotronics  Model 
111A22  with  zero  clamping  power  supply)  located  in  the  nozzle  chamber  behind  the 
nozzle  orifice  (see  figure  2.2).  The  time  resolution  of  the  recorded  pressure  was  1ms 
and  the  pressure  resolution  was  better  than  5  psi. 

Drop  size  measurements  were  recorded  unobtrusively  using  a  single  axis  Aero- 
metric  Phase/Doppler  Particle  Analyzer  (Transmitter  Model  XMT-1100-4S,  Receiver 
Model  RCV-2100,  Motor  Controller  Model  MCB-7100-1,  Signal  Processor  Model 
PDP-3100,  and  Data  Management  System  Model  DMS-4000-5).  This  equipment  ana¬ 
lyzes  the  spatial  and  Doppler  frequency  of  an  interference  pattern  generated  by  two 
intersecting  laser  beams  as  drops  pass  through  the  intersection  point  in  order  to  de¬ 
termine  the  size  and  velocity  of  those  drops  (a  detailed  description  of  this  system  can 
be  found  in  David  Woods’  thesis  [15]).  The  system  was  setup  in  forward  scatter  mode 
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with  the  receiver  oriented  30°  off  the  transmitter  axis.  Drop  size  information  was  re¬ 
corded  continuously  over  many  intermittent  pulses.  That  is,  each  data  run  involves 
numerous  injections  with  the  same  pressure  waveform.  Initiation  of  drop  size  informa¬ 
tion  gathering  was  achieved  through  an  external  trigger  generated  by  the  pressure 
waveform  data  acquisition  system  in  order  to  synchronize  the  two  data  collection  pro¬ 
cesses. 

The  external  trigger  signals  the  start  of  timing  clocks  in  both  the  Keithley  and 
Aerometric  systems.  These  are  the  clocks  which  mark  the  arrival  of  data  with  respect 
to  the  start  of  data  run.  Each  clock  starts  at  time  zero  on  the  trigger  signal,  i.e.  the  time 
of  drop  detection  and  pressure  sampling  are  relative  to  the  initial  trigger  signal. 

Velocities  are  measured  parallel  to  the  axis  of  the  spray  and  are  considered  posi¬ 
tive  travelling  away  from  the  nozzle.  Acceleration  due  to  gravity  is  also  parallel  to  the 
axis  of  the  spray  and  positive  travelling  away  from  the  nozzle. 

2.2  Physical  Properties  of  the  Test  Fluid 

Since  most  studies  have  been  diesel  and  the  primary  application  of  this  work  is 
diesel  injection,  locally  obtained  No.  2  diesel  fuel  was  used  for  these  experiments. 
There  is  no  problems  associated  with  using  diesel  fuel  with  the  hydraulic  equipment; 
however,  it  is  requires  adequate  ventilation. 

The  density  of  the  fuel  was  found  to  be  0.830  g/cm\  Viscosity  was  measured  to  be 
4.464  centistokes  (4.464  x  10 2  kgnv’-s1)  using  a  Cannon  Viscometer  model  OC-B223 
(capillary  tube  efflux  method).  Surface  tension  was  measured  to  be  31.15  dynes/cm 
using  a  DuNouy  Tensiometer  model  70535  (ring  method).  All  tests  were  conducted  at 
the  room  temperature  of  21°c. 

23  Experimental  Procedure 

The  spray  resulting  from  a  given  nozzle  pressure  waveform  varies  from  pulse  to 
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pulse  due  to  the  randomness  of  the  atomization  process.  Averaging  of  the  spray  char¬ 
acteristics  over  a  number  of  spray  events  provides  a  more  general  representation  of 
what  might  be  expected  from  any  given  spray  event.  Thus,  the  connection  between  the 
dynamic  pressure  input  and  the  resulting  spray  is  presented  as  the  average  of  approxi¬ 
mately  seventy  spray  events.  Each  spray  event  is  defined  here  to  be  the  period  start¬ 
ing  10ms  before  the  nozzle  pressure  peak  and  ending  30ms  after  the  pressure  peak.  The 
local  zero  reference  time  from  one  injection  event  to  another  was  picked  to  be  the  time 
of  the  pressure  peak  occurrence  for  each  injection.  Thus,  the  time-resolved  character¬ 
istics  could  be  averaged  by  aligning  each  injection  by  its  peak.  This  assumes  that  each 
waveform  is  nearly  the  same  for  injections  with  peak  pressures  in  the  936  -  1036psi 
range.  The  variation  in  the  pressure  waveform  is  studied  in  section  2.5. 

A  particular  nozzle  pressure  waveform  and  its  spray  ware  analyzed  at  a  particu¬ 
lar  location  by  simultaneously  recording  pressure  and  spray  characteristics  over  a  pe¬ 
riod  of  20  seconds  and  many  injections.  Data  for  individual  pulses  was  then  taken  out 
of  the  continuous  data  stream  with  off-line  processing. 

The  first  priority  was  to  establish  a  continuously  pulsating  spray  with  the  desired 
nozzle  pressure  waveform.  After  the  system  was  pressurized  by  the  supply  pump,  a  si¬ 
nusoidal  signal  was  applied  to  the  control  valve.  The  system  was  allowed  to  run  until 
a  regular  pulsating  spray  was  established  which  takes  a  few  minutes.  It  seemed  that  the 
valve  needs  to  get  warmed  up  before  it  responds  well  to  the  relatively  high  frequency 
applied  to  it.  Once  a  regular  pulsating  spray  was  established,  adjustment  to  the  func¬ 
tion  generator  was  made  while  observing  the  resulting  nozzle  pressure  waveforms.  This 
procedure  was  repeated  for  every  data  run.  Since  there  was  a  separate  data  run  for 
every  new  position,  there  was  the  potential  that  a  single  nozzle  pressure  waveform 
could  not  be  faithfully  repeated  from  one  data  run  to  the  next.  Graphs  presented  in 
section  2.5  show  that  an  acceptable  repeatability  was  maintained  from  one  location  to 
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another  using  the  same  pressure  waveform.  Stopping  the  spray  between  data  runs  was 
done  to  insure  that  the  build  up  of  drops  on  the  chamber  windows  did  not  obstruct  the 
PDPA  system.  The  chamber  windows  were  cleaned  between  data  runs. 

After  establishing  a  pulsating  spray  for  analysis,  approximately  one  thousand  drops 
were  test  recorded  with  the  PDPA  system  to  check  its  operating  parameters.  Typically, 
the  photomultiplier  tube  (PMT)  power  needed  to  be  adjusted  in  order  to  maintain  sen¬ 
sitivity  to  the  smallest  drops  without  incurring  saturation  rejections  [17].  The  size  and 
velocity  range  of  measurement  was  another  important  consideration  from  one  data  run 
to  another.  However,  one  size  and  velocity  range  was  found  to  be  adequate  for  all  of 
the  data  runs. 

Actual  data  runs  are  initiated  quickly  after  adjusting  the  PDPA  system  to  avoid 
excessive  wetting  of  the  chamber  windows.  Starting  the  data  collection  routine  on  the 
Keithley  data  acquisition  system  triggers  simultaneous  recording  of  the  nozzle  pressure 
by  the  Keithley  data  acquisition  system  and  drop  information  by  the  PDPA  system. 
Each  system  records  its  data  independently,  but  each  individual  data  item  in  both  sys¬ 
tems  has  a  time  stamp  relative  to  the  start  of  acquisition.  Thus,  each  data  run  has  a  pair 
of  data  sets.  One  set  contains  the  nozzle  pressure  recorded  every  millisecond  and  the 
other  contains  drop  size  and  velocity  information  recorded  as  the  drops  were  detected. 

Once  the  raw  pressure  and  drop  information  was  collected,  off-line  analysis  was 
then  done  to  remove  data  associated  with  undesired  nozzle  pressure  waveforms  (due 
to  the  supply  pressure  variation)  and  to  perform  the  desired  statistical  functions.  The 
off-line  analysis  references  the  local  pressure  peak  of  each  injection.  Thus,  to  perform 
the  analysis,  the  local  pressure  peaks  must  be  found  in  the  continuous  stream  of  data 
in  a  data  run.  The  number  of  injections  in  a  data  run  dictates  that  the  process  be  au¬ 
tomated.  Local  peaks  are  found  by  using  a  second  order  polynomial  to  interpolate  a 
more  accurate  pressure  peak  from  each  local  peak  in  the  stream  of  pressure  data.  All 
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pressure  pulses  with  interpolated  peaks  in  the  936  to  1036psi  range  were  kept  and  the 
others  rejected  (see  Apendix  A.1).  The  time  offset  between  the  interpolated  peak  and 
the  recorded  peak  pressure  was  also  recorded  in  order  to  accurately  compute  average 
pressure  statistics  for  a  given  data  run  (see  section  2.5,  and  Appendix  A.2). 

Handling  of  the  PDPA  data  required  exporting  the  data  sets  from  Aerometric’s 
binary  format  to  standard  ASCII  format  for  use  with  custom  Fortran  programs.  The 
Aerometric's  software  is  capable  of  handling  a  regular  pulsating  jet;  however,  the  fluc¬ 
tuations  in  pump  supply  pressure  required  editing  of  the  data.  For  that  reason,  the 
Aerometric’s  software  could  not  be  used.  In  addition,  the  use  of  Fortran  programming 
(see  Appendix  A3)  allowed  more  flexibility  in  automating  the  preparation  and  presen¬ 
tation  of  spray  data  in  conjunction  with  the  nozzle  pressure  data. 


Figure  2.3.1  Typical  Gausian  laser  beam  intensity  profile.  Larger 
drops  scaner  more  light  than  small  drops.  Thus,  small  drops  passing 
through  the  periphery  of  the  beam  are  not  sufficiently  illuminated  to 
be  detected  by  the  limited  sensitivity  of  the  receiver. 

There  was  one  problem  associated  with  not  using  the  Aerometric’s  software  for 
analysis.  An  error  in  the  detection  of  small  particles  occurs  when  small  particles  tran- 
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sit  the  edges  of  the  probe  volume.  Since  the  laser  produces  a  beam  with  a  Gausian  in¬ 
tensity  as  shown  in  figure  23.1,  small  particles  will  scatter  less  light  than  larger  particles 
when  they  transit  the  periphery  of  the  probe  volume.  Thus,  the  smallest  size  of  particles 
detected  increases  as  the  transit  path  of  the  particle  moves  away  from  the  center  of  the 
probe  volume.  The  PDPA  software  calculates  the  Gausian  intensity  of  the  probe  vol¬ 
ume  in  order  to  determine  an  actual  probe  volume  for  the  small  drop  sizes.  It  then  cor¬ 
rects  the  distribution  statistics  by  adding  a  number  of  drops  to  the  small  drop  sizes  in 
order  to  compensate  for  the  different  actual  probe  volumes. 

The  raw  data  exported  from  the  PDPA  software  does  not  contain  corrections  for 
the  different  actual  probe  volumes,  but  rather  only  gives  the  data  on  the  actual  drops 
detected.  In  order  to  correct  my  ‘custom’  distributions,  the  correction  factors  were  com¬ 
puted  from  each  size  range  in  the  PDPA  generated  distributions  (see  Appendix  A.4). 
These  correction  factors  were  then  used  to  correct  the  custom  distributions  for  varia¬ 
tion  in  actual  probe  volumes  by  adding  a  proportional  number  of  drops  throughout  the 
range  of  my  ‘custom’  distributions. 

2.4  Sample  Locations 

Figure  2.4.1  shows  the  six  fixed  drop  size  and  velocity  sampling  locations  used. 
Sampling  in  these  locations  is  achieved  by  moving  the  nozzle  relative  to  the  probe  vol¬ 
ume.  Since  the  spray  chamber  does  not  move  with  the  nozzle,  there  is  the  potential  that 
the  change  from  one  sampling  location  to  another  will  influence  the  spray  dynamics. 
Observation  of  the  spray  in  each  location  showed  no  visible  deflection  of  the  spray 
caused  by  the  chamber  floor  or  walls. 

2.5  Study  of  the  Nozzle  Pressure  Waveforms  Used 

The  acceleration  of  the  liquid  in  the  atomizing  jet  is  related  to  the  nozzle  pres¬ 
sure  rate  of  change.  Several  smooth  pressure  waveforms  were  used  to  create  pulsat¬ 
ing  jets.  This  allowed  the  observation  of  the  change  in  spray  characteristics  of  the  liq- 
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uid  jet  undergoing  different  acceleration  patterns.  Both  symmetric  and  asymmetric 
pressure  waveforms  were  created.  Variation  of  the  period  of  the  symmetric  waveforms 
produced  approximately  equal  changes  in  the  leading  and  trailing  edges  of  the  time- 
resolved  pressure  waveform.  The  system  was  also  able  to  significantly  vary  the  trail¬ 
ing  edge  of  the  waveform  yielding  asymmetric  pressure  pulses.  However,  the  system 
did  not  accommodate  changing  the  leading  edge  of  the  pressure  waveform  except  in 
the  symmetric  case. 

Since  the  spray  data  presented  is  based  upon  the  average  of  many  pulses,  it  is  ap¬ 
propriate  to  consider  the  variation  from  one  pulse  to  another.  Recall  that  the  pressure 
peak  was  determined  by  second  order  polynomial  interpolation  and  that  there  is  a  time 
offset  of  the  interpolated  peak  from  the  measured  peak.  For  averaging,  each  pulse  is 
time-aligned  by  its  interpolated  peak.  Then  pressure  values  are  linearly  interpolated 
at  the  peak  offset  every  millisecond  from  10ms  before  the  peak  to  30ms  after  the  peak. 

The  average  nozzle  pressure  waveforms  for  each  data  run  are  presented  in  sets. 
Comparison  of  the  same  pressure  waveform  reflects  the  need  to  verify  the  pressure 
waveform  for  each  data  run  after  shutting  down  the  pulsating  spray  system  between 

Nozzle  EM  X  (radial  direction) 

A*' 

J ►  . 
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A-B  (0.32cm  ,8.4cm)  •  ;  *B-B  (1.32cm,  8.4cm) 


B-A  (1.32cm, 23.9cm) 

A-A  (0.32cm, 23.9cm)  •••  *D-A  (232cm,  23.9cm) 

C-A  (0.82cm, 23.9cm) 

*  Z  (axial  direction) 

Figure  2.4.1  Spray  sampling  locations  downstream  of  the  nozzle. 
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runs.  Table  5.3.1  gives  the  results  of  calculating  the  standard  deviation  for  three  typi¬ 
cal  data  runs. 


Results  of  the  spray  analysis  will  refer  to  particular  pressure  waveforms  as  they 
are  presented  here  in  the  following  graphs,  figure  2.5.1.  There  are  three  symmetric 
pressure  waveforms  denoted  by  their  approximate  duration  of  6, 14,  and  16ms.  In  ad¬ 
dition,  there  are  two  asymmetric  pressure  waveforms  denoted  gentle  trailing  edge 
(Gentle  TE)  and  steep  trailing  edge  (Steep  TE). 


Pulse 

Waveform 

Average 

Standard 

Deviation* 

Maximum 

Standard 

Deviation 

Standard  Deviation  of 
Pressure  Peak 

I6mt 

Symmetric  at 
&*A 

±1B6p» 

±49  Ipei 

±!6.2pei 

Gentle  TE  at 
B-A 

±24  4pei 

±65.8p« 

±9  Ipei 

Steep  TE  at 
B-A 

±2V6psi 

±69.7psi 

±16.4p» 

Table  2.5.1  Typical  nozzle  pressure  waveform  statistics. 
*  The  average  standard  deviation  is  the  average  of  standard 
deviations  computed  on  the  pressure  values  at  each 
millisecond. 
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TIME  (ms.  RELATIVE  TO  THE  PRESSURE  PEAK) 


6,14,  and  16ms  symmetric  pulses  at  location  A*B  Gentle  TE  pulses  at  I  motions  A-B,  B-B,  B-A, 

and  A-A 


16ms  symmetric  pulses  at  locations  A-B,  B-B, 
B-A,  and  A-A 


Steep  TE  pulses  at  locations  A-B,  B-B,  B-A, 
and  A-A 


16ms  symmetric  pulses  at  locations  A-A,  C-A,  B-  Comparison  of  16ms  symmetric,  gentle  TE,  and 
A,  and  D-A  steep  TE  pulses  at  location  B-A 


Figure  2.5.1  Nozzle  pressure  waveforms  used  for  data  runs  at  particular  locations.  Pressure 

represents  the  average  pressure  of  all  the  spray  pulses  used  in  a  data  run. 
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Chapter  3 


Results 


3.0  General 

Drop  size  and  velocity  information  is  presented  as  the  average  of  measurements 
taken  over  approximately  70  spray  events.  In  similar  work,  Koo  and  Martin  [5]  found 
consistent  results  from  data  run  to  data  run  using  50  to  100  spray  events.  Thus,  where 
data  is  presented  over  the  course  of  a  spray  event,  it  is  actually  the  average  of  data  re¬ 
corded  for  a  number  of  spray  events.  The  exact  number  of  spray  events  used  depends 
upon  the  number  of  pulses  accepted  (see  section  2.3  for  details  of  the  acceptance  cri¬ 
teria)  in  a  given  data  run. 

Comparison  of  results  from  one  pressure  waveform  to  another,  such  as  the  6ms 
symmetric  to  the  16ms  symmetric,  requires  discretion.  There  is  a  difference  in  the  du¬ 
ration  of  the  spray  event  from  one  waveform  to  another.  This  will  influence  the  over¬ 
all  number  of  drops  per  pulse  because  of  the  change  in  injected  quantity.  However,  the 
change  in  efficiency  of  atomization  is  more  important  to  observe,  and  this  change  can 
be  seen  in  the  change  of  distribution  profiles. 

3.1  Time  Integrated  Distribution  of  Drop  Size  and  Velocity 

Characteristics  of  the  spray  produced  by  various  pressure  waveforms  in  the  nozzle 
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were  integrated  over  the  duration  of  the  spray  event.  The  spatial  location  of  the  sam¬ 
pling  is  fixe*1  during  the  course  of  the  integrations.  This  data  is  the  average  total  drop 
distribution  falling  in  a  specific  location  relative  to  the  nozzle. 

Drop  sizes  of  all  of  the  drops  recorded  in  the  sample  volume  during  each  spray 
event  are  tabulated  and  grouped  (by  size)  into  one  micron  groups  from  1  to  75  microns. 
The  totals  are  then  divided  by  the  number  of  selected  spray  events,  yielding  the  drops 
per  pulse  for  comparison,  and  the  results  are  graphed  in  figure  3.1.1. 

None  of  the  sprays  have  a  monodispersed  drop  size  distribution.  Regaidless  of 
pressure  waveform  and  sampling  location,  each  total  size  distribribution  has  its  peak 
drops  per  pulse  occuring  in  the  10-20  micron  size  range.  The  number  of  drops  per  pulse 
falls  off  sharply  in  the  0-10  micron  range  ar.d  gradually  from  20  to  75  microns. 

Comparing  the  three  symmetric  waveforms  at  location  B-A,  we  see  an  approxi¬ 
mately  equal  rise  in  drops  per  pulse  from  0  to  10  microns.  However,  from  10  to  20  mi¬ 
crons  there  is  some  irregular  separation  in  the  distributions,  and  from  20  to  75  microns 
each  distribution  declines  at  approximately  the  same  rate  though  separate.  The  pres¬ 
sure  waveform  with  the  longest  duration  (16ms)  has  the  consistently  highest  drops  per 
pulse  while  the  shortest  (6ms)  has  the  lowest  drops  per  pulse 

One  symmetric  waveform  (16ms)  was  also  studied  in  five  spatial  locations  (A-B, 
B-B,  A-A,  C-A,  and  D-A)  in  addition  to  B-A.  Two  upstream  locations,  A-B  and  B-B, 
are  compared  with  two  parallel  downstream  locations,  A-A  and  B-A.  We  can  see  that 
the  spray  is  much  more  dense  at  the  near  axial  location  A-B  by  its  broad,  high  drops 
per  pulse  distribution.  Location  B-B  shows  a  low  level  throughout  the  range,  which  is 
probably  due  to  the  spray  cone  not  reaching  out  that  far.  The  downstream  locations 
have  a  generally  lower  level  of  drops  than  location  A-B  throughout  the  size  range  due 
to  the  spread  in  the  spray  cone.  In  addition,  the  rise  in  drops  per  pulse  from  1  to  10 
microns  is  less  sharp  downsteam  than  at  the  location  A-B,  upstream. 
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Figure  3.1.1  Time-integrated  drop  size  distributions  for  different  pressure  waveforms  and  sampling 
locations. 


The  16ms  Symmetric  pressure  pulse  was  also  compared  in  the  four  radial  down¬ 
stream  locations,  A-A,  C-A,  B-A,  and  D-A.  Outer  locations  B-A  and  D-A  show  that 
there  are  more  small  drops  per  pulse  transiting  these  locations  than  the  inner  locations. 

The  Gentle  TE  pressure  waveform  was  compared  at  the  four  spatial  locations  A- 
B,  B-B,  B-A,  and  A-A.  Again,  the  highest  drops  per  pulse  throughout  the  size  range 
were  found  at  location  A-B  although  the  distribution  has  a  more  dramatic  peak  in  the 
10-20  micron  range  than  the  16ms  Symmetric  pulse.  Downstream  locations  have  a  less 
steep  rise  in  drops  per  pulse  from  1  to  10  microns  than  A-B.  They  are  also  at  higher 
levels  than  the  16ms  Symmetric  pulses  at  the  same  locations.  The  outer  location  B-A 
has  more  drops  per  pulse  in  the  10-20  micron  range  than  the  near  axial  location  A-A 
as  we  see  with  the  16ms  Symmetric  pulse. 

The  Steep  TE  pressure  waveform  was  compared  at  the  four  spatial  locations  A- 
B,  B-B,  B-A,  and  A-A.  We  see  that  the  upstream  location  A-B  has  a  relatively  low 
drops  per  pulse  throughout  the  size  range  in  contrast  with  the  16ms  Symmetric  and 
Gentle  TE  pulses.  The  downstream  outer  location  B-A  has  broad  peak  from  5  to  20 
microns  and  has  a  higher  drops  per  pulse  throughout  the  whole  size  range  than  the 
near  radial  location  A-A. 

Note  the  consistent  spike  in  the  number  of  drops  per  pulse  at  43  microns,  which 
is  typically  25%  greater  than  adjacent  drops  per  pulse  levels.  Other  spikes  or  fluctua¬ 
tions  occur,  but  they  do  not  appear  consistently  and  are  probably  due  to  statistical  con¬ 
siderations.  The  repeated  occurrance  of  this  spike  at  43  microns,  without  regard  to 
pressure  waveform  and  sampling  location,  indicates  that  it  is  probably  a  physical  phe¬ 
nomena  and  not  an  error.  Testing  of  other  sprays  showed  no  increased  sensitivity  of 
the  Aerometric’s  system  to  this  size  drop. 

Velocities  of  all  of  the  drops  recorded  in  the  sample  volume  during  each  spray 
event  are  also  tabulated  and  grouped  into  lm/s  groups  from  1  to  75m/s.  The  resulting 
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Figure  3.1.2  Time-integrated  drop  velocity  distributions  for  different  pressure  waveforms  and  sampling 
locations. 


21 


totals  are  then  divided  by  the  number  of  spray  events,  yielding  the  drops  per  pulse  in 
figure  3.1.2  for  comparison.  There  is  quite  a  bit  more  fluctuation  in  the  individual  time- 
integrated  velocity  profiles. 

The  velocity  profiles  are  bimodal  at  all  downstream  locations  for  every  pressure 
waveform.  Peak  drops  per  pulse  occur  at  approximately  10-15  m/s  and  25-30  m/s.  The 
10-15  m/s  peak  surpasses  the  25-30  m/s  peak  for  all  the  symmetric  pressure  waveforms 
and  the  Steep  TE  waveform,  but  the  opposite  is  true  for  the  Gentle  TE  waveform. 

Upstream,  the  velocity  profiles  for  the  16ms  Symmetric  and  Gentle  TE  pressures 
are  broad  near  the  axis  at  A-B  and  sharp  away  from  the  axis  at  B-B.  The  peak  drops 
per  pulse  are  in  the  same  magnitude  range  as  the  downstream  locations  except  for  the 
Steep  TE  pressure  waveform  which  has  much  lower  peaks  in  contrast. 

3.2  Temporally  and  Spatially  Integrated  Drop  Size  Distribution 

Drop  size  information  was  integrated  over  the  course  of  the  spray  event  at  radial 
locations  A-A,  C-A,  B-A,  and  D-A  for  the  16ms  Symmetric  pressure  waveform.  The 
total  drops  are  divided  by  the  number  of  spray  events. 
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Here  we  get  a  more  refined  view  of  the  total  downstream  drop  distribution  than 
looking  at  one  location  only.  Common  characteristics  of  the  respective  individual  pro¬ 
files  in  figure  3.1.1  can  be  found  in  the  total  combination  of  the  graphs.  For  instance, 
there  is  the  steep  rise  in  the  drops  per  pulse  from  0  to  10  microns.  The  peak  drops  per 
pulse  occurs  in  the  10-20  micron  range,  and  the  drops  per  pulse  falls  off  gently  from 
20  to  75  microns.  The  spike  at  43  microns  is  also  very  clear. 

33  Time  Evolution  of  Drop  Size  and  Velocity  Distributions 

The  size  of  all  of  the  drops  recorded  in  the  sample  volume  during  each  2ms  pe¬ 
riod  in  the  spray  event  are  tabulated  in  one  micron  groups  from  0-75  microns;  while 
the  velocity  of  all  the  drops  recorded  are  tabulated  in  lm/s  groups  from  0-75  m/s.  The 
resulting  totals  are  then  divided  by  the  number  of  spray  events  in  order  to  get  an  av¬ 
erage  number  of  drops  per  pulse  for  comparison.  Each  set  of  totals  produces  an  “in¬ 
stantaneous'  distribution  histogram  for  every  other  millisecond  of  the  spray  event. 

There  are  a  total  of  fifty-eight  drop  size  and  velocity  histograms  to  cover  a  spray 
event  with  2ms  resolution.  Therefore,  only  the  distribution  histograms  for  the  16ms 
symmetric  pressure  pulse  at  location  B-A  are  presented  here  (see  figures  3.3. 1,2). 

Both  the  size  and  velocity  histograms  show  that  the  spray  does  not  arrive  at  the 
downstream  location  B-A  until  the  2ms  period  after  the  pressure  pulse  reaches  its 
peak.  Drops  continue  passing  this  location  46ms  after  the  pressure  peak. 

The  size  histograms  start  out  broad  when  the  spray  reaches  the  sampling  location, 
but  the  small  drops  in  the  10-20  micron  range  dominate  quickly  until  about  40ms  af¬ 
ter  the  pressure  peak.  There  is  a  small  bimodal  indication  which  is  most  pronounced 
12-14ms  after  the  pressure  peak.  The  bimodal  peaks  are  at  roughly  16  and  40  microns. 
Although  it  is  impossible  to  determine  the  exact  correspondance  of  the  nozzle  pres¬ 
sure  with  the  drops  found  downstream  of  the  nozzle,  the  bimodal  indication  occurs 
during  the  time  when  the  peak  number  of  drops  is  being  recorded. 
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Axial  velocity  of  the  drops  also  starts  out  broad,  but  becomes  quite  narrow  in 
range.  The  velocity  distributions  from  4  to  12ras  after  the  pressure  peak  have  nearly 
symmetric  distributions  around  30ra/s.  After  12ms  the  velocity  distributions  tend  to¬ 
wards  a  narrow  band  at  lOm/s. 

Average  data  is  extracted  from  distribution  histograms  for  each  waveform  at  each 
location  to  yield  graphs  of  Sauter  mean  diameter,  drops  per  pulse,  and  mean  velocity 
(see  figures  33.3-5).  These  graphs  summarize  the  instantaneous  distributions  histo¬ 
grams.  The  time  resolution  of  these  summary  graphs  is  1ms  compared  to  the  two  mil¬ 
lisecond  resolution  presented  in  figs.  33.1  and  33.2.  Original  histogram  data  has  1ms 
resolution,  but  histograms  were  combined  to  reduce  the  number  of  histograms  pre¬ 
sented  in  half. 

3.4  Continuous  Jets 

Drop  size  and  velocity  data  were  recorded  for  continuous  jets  of  789, 888,  and  986 
psi  at  location  B-A.  Since  these  are  continuous  jets,  the  convenient  measure  of  drops 
per  pulse  cannot  be  used.  Rather,  the  number  of  drop  totals  are  presented.  Compari¬ 
son  can  be  made  to  the  distributions  for  pulses  by  qualitative  evaluation  of  the  distri¬ 
bution  profiles.  The  size  and  velocity  mean  data  (see  table  3.4.1)  can  be  directly  com¬ 
pared  to  those  means  presented  for  pulsating  sprays. 


Table  3.4.1  Statistics  for  continuous  sprays  at  789, 888,  and 
986psi  nozzle  pressures. 
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In  figure  3.4.1  we  see  at  888  and  986  psi  that  the  drop  size  distribution  rises  steeply 
from  0  to  10  microns,  peaks  in  the  10-20  micron  range,  and  declines  gently  from  20  to 
75  microns.  All  three  pressure  ranges  have  a  spike  at  43  microns.  Other  significant 
spikes  occur  at  7  and  23  microns  for  the  888psi  profile  and  12  microns  for  the  986  psi 
profile.  The  789  psi  drop  size  distibution  is  broad  with  a  gentle  peak  in  the  30-45  mi¬ 
cron  range  except  for  a  spike  at  43  microns. 

The  velocity  profiles  are  somewhat  symmetric.  We  can  see  that  the  profiles  for  888 
and  986  psi  are  quite  close  with  peaks  in  the  27-33  m/s  range.  The  789  psi  velocity  pro¬ 
file  peak  is  in  the  22-30  m/s  range. 
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Figure  3 33  Time  evolution  of  the  Sauter  mean  diameter  for  different  pressure  waveforms  and 
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Figure  3 3A  Time  evolution  of  the  average  drops/pulse  recorded  for  different  pressure  waveforms 
and  different  sampling  locations. 
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Figure  3 35  Time  evolution  of  the  velocity  recorded  for  different  pressure  waveforms  and  different 
sampling  locations. 
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Chapter  4 


Discussion 


4.1  Observed  Trends 

Many  trends  in  the  spray  characteristics  produced  by  different  nozzle  pressure 
waveforms  are  identified  from  the  results  presented  in  Chapter  3.  These  trends  show 
the  influence  of  a  dynamic  pressure  input  upon  the  drop  size  and  velocity  distributions, 
as  well  as  the  propagation  of  the  drops  from  the  nozzle. 

Comparing  the  time-integrated  drop  size  distributions  at  location  B-A  (see  figure 
3.1.1),  three  drop  size  ranges  are  considered.  The  production  of  0  to  10  micron  drops 
is  not  dependent  upon  the  pressure  waveform  since  all  of  the  distributions  are  nearly 
the  same  in  this  range.  From  10  to  20  microns,  where  each  distribution  peaks,  drop 
production  is  quite  dependent  upon  the  pressure  waveform.  In  the  largest  size  drops, 
from  20  to  75  microns,  the  differential  in  drop  production  between  the  different  pres¬ 
sure  waveforms  is  consistent  with  injection  duration.  That  is,  the  longer  the  duration 
of  the  pulse  the  more  large  size  drops  are  produced. 

Location  B-A  was  used  to  make  comparisons  between  the  various  nozzle  pressure 
waveforms.  All  of  the  time-integrated  size  distributions  at  this  location  show  that  drops 
in  the  10  to  20  micron  diameter  range  were  produced  the  most.  The  common  feature 
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of  each  pressure  waveform  was  a  peak  pressure  of  approximately  986  psi.  Thus,  it  is 
generalized  that  this  particular  nozzle  will  produce  mainly  drops  in  the  10  to  20  micron 
diameter  range  when  the  peak  pressure  reaches  approximately  986  psi  (maintaining 
identical  fluid  and  temperature)  for  any  pressure  waveform. 

Among  symmetric  pulses,  the  connection  between  peak  pressure  and  production 
of  drops  in  the  10  to  20  micron  size  range  is  reinforced.  Since  the  symmetric  pulses 
were  nearly  sinusoidal,  the  period  of  time  the  pressure  remains  near  peak  is  deter¬ 
mined  by  the  period  of  the  pulse.  Longer  symmetric  pulses  have  a  longer  period  near 
peak  pressure;  thus,  these  pulses  should  produce  more  drops  in  the  10  to  20  micron 
range.  In  fact,  the  symmetric  waveform  with  the  longest  duration,  the  16  ms  pulse,  pro¬ 
duced  the  most  drops  within  this  range.  Likewise,  the  opposite  is  also  observed.  The 
shortest  pulse  (6  ms)  produced  the  least  drops  in  this  range. 

Comparing  the  symmetric  pulse  with  the  two  asymmetric  pulses  also  supports  the 
link  between  peak  pressure  and  production  of  drops  in  the  10  to  20  micron  size  range. 
The  Gentle  TE  pulse  clearly  has  the  longest  period  near  peak  pressure  (see  figure 
2.5.1),  and  it  has  the  highest  production  of  drops  in  this  range.  The  distribution  for  the 
Steep  TE  pulse  shows  more  drops  in  this  size  range  than  the  14  ms  Symmetric  pulse, 
but  neither  has  an  observable  difference  in  pulse  shape  near  peak  pressure. 

A  link  is  made  between  the  production  of  larger  drops  and  the  trailing  pressure 
function  in  the  nozzle.  The  sharp  decline  in  nozzle  pressure  after  the  peak  reduces  the 
number  large  drops  found  at  downstream  location  B-A.  Comparing  the  14  ms  Symmet¬ 
ric  with  the  two  asymmetric  pulses,  the  Steep  TE  pulse  has  the  fewest  drops  in  the  larg¬ 
est  size  range.  The  sharp  pressure  decline  in  the  Steep  TE  pulse  reduces  the  period  of 
injection  where  larger  drop  production  is  established. 

It  may  also  be  true  that  a  sharp  rise  in  pressure  to  the  peak  also  reduces  the  num¬ 
ber  of  large  drops  found  at  downstream  location  B-A.  Comparing  the  symmetric  pulses, 
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the  6  ms  pulse  has  the  fewest  drops  in  the  larger  sizes,  but  this  does  not  necessarily 
prove  that  a  sharp  decline  in  trailing  nozzle  pressure  produces  less  large  drops  as  stated 
in  the  previous  paragraph.  The  leading  pressure  function  matches  the  trailing  pressure 
function  with  symmetric  pulses.  Thus,  a  sharp  rise  in  pressure  may  also  mean  less  large 
drops  are  produced.  A  sharp  rise  in  pressure  certainly  would  reduce  the  period  where 
production  of  large  drops  could  take  place  at  the  onset  of  injection,  and  may  inhibit 
the  production. 

A  one  and  two-dimensional  spatial  comparison  is  made  from  the  time-integrated 
drop  size  distributions.  The  four  primary  locations  A-A,  A-B,  B-B,  B-A  are  intended 
to  give  a  global  indication  of  spray  activity.  The  two  additional  locations  C-A  and  D- 
A  for  the  16  ms  Symmetric  pulse  in  combination  with  A-A  and  B-A  give  more  detailed 
insight  into  spray  activity  along  a  radial  line  from  the  axis  of  the  spray. 

Location  A-B  is  closest  to  the  nozzle  exit.  The  time-integrated  size  distributions 
vary  with  the  input  pressure  waveform.  Both  the  16  ms  Symmetric  and  Gentle  TE 
pulses,  which  have  relatively  long  durations,  produce  their  largest  distributions  at  this 
location.  There  are  two  size  ranges  present  in  these  distributions,  0-22  microns  and  22- 
45  microns.  The  smaller  p»nge  is  the  peak  range;  while  the  larger  range  is  a  broad  pla¬ 
teau. 

The  Steep  TE  pulse  produced  a  relatively  low  distribution  at  A-B  compared  with 
its  distributions  dow  nstream.  The  sharp  cutoff  in  pressure  of  the  Steep  TE  waveform 
does  not  allow  the  liquid  jet  to  breakup  as  effectively  until  further  downstream. 

Very  few  drops  are  found  at  location  B-B  for  any  pressure  waveform.  This  loca¬ 
tion  is  radial  to  A-B,  and  the  lack  of  drops  indicates  that  the  spray  angle  is  less  than 
9®.  However,  a  few  drops  of  all  sizes  are  found  there,  probably  due  to  recirculation  or 
sputtering  from  the  nozzle  face. 
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Downstream,  at  locations  A-A  and  B-A,  the  distribution  profiles  are  similar  for 
all  pressure  waveforms.  However,  more  drops  are  found  in  the  outer  location  B-A.  A 
comparison  of  radial  locations  is  made  for  the  16  ms  Symmetric  pulse,  which  helps 
clarify  this  observation. 

We  see  from  the  radial  comparison  of  sampling  locations  a  step-like  difference 
between  the  time-integrated  size  distribution  at  C-A  and  the  one  at  B-A.  A  higher 
number  of  drops  in  the  peak  range  for  the  outer  locations  is  the  difference.  That  is, 
there  are  more  small  drops  in  the  outer  locations  than  in  the  inner  two  locations.  It  is 
difficult  to  conclude  why  we  see  this  trend.  Since  figure  3.4.4  shows  that  the  bulk  of  the 
drops  arriving  at  outer  location  B-A  before  the  inner  location  A-A,  more  of  the  drops 
at  B-A  must  be  from  the  increasing  nozzle  pressure.  The  spray  cone  must  be  expand¬ 
ing  to  the  point  where  more  drops  are  found  in  the  periphery  than  along  the  axis  by 
the  time  the  spray  reaches  the  downstream  distance  of  locations  A-A  and  B-A. 

These  radial  locations  are  also  used  to  give  an  idea  of  the  total  drop  size  distri¬ 
bution  of  the  spray.  Figure  3.2.1  totals  the  drops  from  locations  A-A,  C-A,  B-A,  and 
D-A  for  a  coarse  representation  of  the  distribution  which  might  be  expected  from  a 
complete  cross-sectional  sampling.  The  distribution  closely  resembles  the  distribution 
for  the  constant  888  psi  spray  (see  figure  3.4.1). 

The  time-integrated  size  distributions  are  supplemented  by  time-integrated  veloc¬ 
ity  distributions  (see  figure  3.1.2).  As  reported  in  the  results  section,  these  velocity  dis¬ 
tributions  were  bimodal  at  all  downstream  locations,  broad  at  location  A-B,  and  sharply 
peaked  at  location  B-B. 

The  bimodal  velocity  distributions  taken  alone  are  hard  to  decipher.  However, 
they  point  out  that  something  is  taking  place  between  upstream  location  A-B  and  the 
downstream  locations,  since  the  distribution  is  not  bimodal  at  A-B.  One  possible  ex¬ 
planation  is  that  a  denser  sheet  of  drops  moves  out  radially  downstream.  It  then  col- 
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lapses  as  the  pressure  declines.  The  sheet  passes  the  sampling  locations  downstream 
twice  with  different  average  velocities  each  time. 

Another  explanation  of  the  biraodal  velocity  distributions  downstream  is  that  coa¬ 
lescence  produces  a  peak.  Smaller  drops  pass  through  the  sampling  region  with  a  high 
average  velocity  producing  one  peak.  However,  significant  coalescence  may  also  pro¬ 
duce  a  sub-group  of  drops  with  a  lower  average  velocity  resulting  in  another  peak. 

The  time-integrated  velocity  distribution  at  B-B  is  not  bimodal  for  any  pressure 
waveform.  There  is  one  sharp  peak  at  approximately  8  m/s.  Thus,  these  drops  are  prob¬ 
ably  recirculating  or  freefall  as  previously  mentioned. 

In  contrast  to  the  time-integrated  distributions,  the  time-resolved  distributions 
give  an  evolutionary  picture  of  the  injection  event.  Since  there  is  some  distance  from 
the  nozzle  pressure  sampling  location  and  the  drop  sampling  locations,  there  is  no  di¬ 
rect  link  between  instantaneous  pressure  and  instantaneous  drop  distribution  charac¬ 
teristics.  However,  we  can  make  some  observations  of  the  instantaneous  drop  distri¬ 
bution  characteristics,  and  draw  some  conclusions. 

There  is  a  pattern  consistent  through  all  of  the  Sauter  mean  diameter  (d32)  graphs 
(see  figure  3.3.3).  The  d32  generally  stays  within  a  range  of  0  to  50  microns.  Dramatic 
fluctuations  of  d32  occur  during  periods  where  there  are  few  drops.  These  fluctuations 
probably  represent  both  the  inconsistent  drop  size  production  and  statistical  inadequa¬ 
cies  of  low  drop  numbers.  Comparison  of  the  d32  graphs  with  their  respective  drops 
per  pulse  graph  shows  that  the  d32  is  most  stable  during  the  period  of  high  drop  per 
pulse  rates. 

A  distinct  minimum  d32  occurs  during  the  stable  period.  Since  there  is  a  time  lag 
between  the  pressure  waveform  and  the  events  downstream  in  the  spray,  it  is  difficult 
to  link  the  two  observations.  However,  this  minimum  d32  occurs  with  the  peak  aver- 
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age  velocity.  Since  the  average  velocity  follows  the  input  pressure,  we  conclude  that  the 
d32  has  an  inverse  relationship  with  pressure. 

The  Sauter  mean  diameter  at  location  B-A  downstream  stays  within  35-45  microns 
during  the  stable  period  for  each  of  the  tested  pressure  waveforms.  Upstream,  near  the 
spray  axis  (A-B),  d32  was  smaller  than  either  of  the  sample  locations  downstream. 
Radial  sample  locations  downstream  showed  no  significant  variation  in  the  d32  pat¬ 
tern.  Smaller  drops  at  upstream  location  A-B  rather  than  downstream  at  A-A  and  B- 
A  indicates  that  coalescence  is  taking  place  axially. 

The  time-evolution  of  the  drops  per  pulse  is  an  indicator  of  spray  activity  and/or 
density  in  the  sample  region.  In  general  the  drops  per  pulse  level  rises  very  fast,  lev¬ 
els  out,  then  falls  quickly  for  a  short  time,  and  finally  falls  slowly  towards  zero.  The 
longer  the  pulse,  the  longer  the  period  where  the  drops  per  pulse  remain  at  a  high 
level. 

Significant  numbers  of  drops  arrive  downstream  in  the  near  axial  sampling  loca¬ 
tion  A-A  later  than  in  the  outer  location  B-A.  Then  20  ms  after  the  pressure  peak  the 
number  of  drops  occurring  in  radial  positions  A-A,  C-A,  B-A,  and  D-A  appears  equal. 
Note  that  drop  sizes  are  at  their  lowest  at  this  time  and  velocity  has  leveled  off  to  9  m/ 
s.  Thus,  these  are  probably  free-fall  drops  that  are  recirculating,  since  they  are  smaller 
than  Rayleigh  mode  drops. 

Earlier  we  discussed  the  bimodal  time-integrated  velocity  distributions  at  loca¬ 
tions  A-A  and  B-A.  When  we  look  at  the  time-resolved  evolution  of  average  drop  ve¬ 
locity  in  figure  3.3.5,  we  see  bimodal  characteristics  at  the  upstream  location  A-B  and 
downstream  location  A-A.  This  supports  the  previous  idea  that  the  propagation  of  the 
spray  is  in  the  form  of  a  dense  sheet  moving  out  ward  radially  and  then  collapsing. 
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4.2  Correlation  to  Continuous  Spray 

Using  the  previous  assumption  that  the  minimum  in  the  evolution  of  the  Sauter 
mean  diameter  for  intermittent  sprays  results  from  the  maximum  nozzle  pressure,  we 
see  that  Sauter  mean  diameter  of  the  continuous  spray  at  888  and  986  psi  (see  figures 

3.4.1)  match  Sauter  mean  diameters  in  the  intermittent  sprays  (see  figure  33 3)  at  simi¬ 
lar  pressures.  Moreover,  the  corresponding  average  drop  velocities  also  match. 

Note  that  the  Sauter  mean  diameter  dramatically  increases  and  the  mean  veloc¬ 
ity  decreases  in  the  nozzle  pressure  drop  from  888  to  789  psi  in  the  continuous  jets.  Yet, 
there  is  no  is  no  appreciable  change  in  the  spray  characteristics  between  the  888  and 
986  psi. 

Looking  at  the  intermittent  sprays  for  a  similar  change  as  the  pressure  increases 
or  decreases  from  approximately  900  to  800  psi  is  difficult  because  no  direct  time  phase 
relationship  is  established  between  pressure  sampling  and  drop  sampling.  However,  in 
general  we  do  not  see  the  dramatic  change  in  size  during  the  evolution  of  the  intermit¬ 
tent  sprays.  There  does  seem  to  be  a  characteristic  steep  decline  in  average  velocity 
following  the  peak  plateau  in  the  evolution  of  average  velocities  for  intermittent  sprays. 

43  Special  Observation 

There  is  a  characteristic  peak  of  drops  produced  at  43  microns  which  shows  up  in 
the  time-integrated  drop  size  distributions  of  the  intermittent  sprays  (see  figures  3.1.1 
and  3.2.1).  This  peak  is  very  narrow  in  size  range  and  is  typically  25%  greater  than  ad¬ 
jacent  size  ranges.  It  is  clear  that  the  nozzle  produces  this  drop  size  with  greater  effi¬ 
ciency  than  others.  Since  the  spike  occurs  with  every  pressure  waveform  there  must  be 
a  common  condition  that  occurs  with  each  waveform,  or  the  spike  is  not  waveform  de¬ 
pendent.  Since  see  the  spike  in  the  continuous  pressure  size  distributions  (see  fig. 

3.4.1) ,  we  can  conclude  that  the  spike  is  not  a  function  of  a  dynamic  pressure  input. 

Interestingly,  this  spike  is  very  close  the  Sauter  Mean  Diameter  of  the  888  and  986 
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psi  continuous  sprays.  It  seems  that  the  nozzle  has  an  optimum  efficiency  at  45  microns 
for  continuous  sprays.  The  dynamic  pressure  inputs,  although  reaching  a  peak  pressures 
similar  to  the  continuous  sprays,  shift  the  efficiency  to  the  smaller  10-20  micron  range. 

Going  back  to  theory  [16],  the  growth  of  many  disturbances  of  different  frequen¬ 
cies  is  responsible  for  atomization  with  one  frequency  dominating.  In  the  identifica¬ 
tion  of  the  spike  at  43  microns,  it  was  also  matKtioned  that  there  were  other  spikes 
present.  The  dynamic  pressure  input  must  enhance  s  competing  frequency  in  order  to 
overwhelm  the  geometrically  advaMaged  bt.qv.enay  producing  43  micron  drops,  since 
it  produces  10-20  micron  drops  primarily. 
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Chapter  5 


Conclusion 


In  the  field  of  intermittent  liquid  jets,  there  is  a  definite  influence  on  the  drop  sizes 
and  velocities  by  a  dynamic  pressure  input  to  the  nozzle.  This  study  shows  that  by 
changing  the  ascent  and  decline  of  the  nozzle  pressure,  while  maintaining  the  same 
peak  pressure,  different  results  are  achieved.  Moreover,  these  results  vary  from  just  a 
simple  change  in  injection  volume.  There  are  changes  in  the  spray  drop  size  and  ve¬ 
locity  distributions. 

In  addition,  evidence  has  been  presented  which  indicates  the  nozzle  used  can  pro¬ 
duce  one  drop  size  with  greater  efficiency  than  other  sizes  at  high  injection  pressures. 
In  particular,  the  appearance  of  the  43  micron  drop  size  spike  warrants  investigation. 
Further  work  may  validate  theories  predicting  a  disturbance  of  a  given  frequency  that 
will  dominate  atomization  under  the  conditions  used.  To  do  continue  on  this  line  of 
study  requires  the  measurement  of  instantaneous  flow  rate  which  was  unavailable  for 
this  work.  Instantaneous  flow  rate  would  give  instantaneous  average  velocity  across  the 
nozzle  exit;  thus,  a  Reynold’s  number  could  be  computed. 

The  broad  nature  of  this  work  lends  itself  very  well  to  more  detailed  studies. 
There  are  two  directions  in  which  to  explore.  The  first  direction  is  to  study  the  influ- 
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ence  of  the  nozzle  pressure  waveform  with  more  variation.  The  second  is  to  explore 
the  spray  produced  with  a  given  waveform  with  more  resolution.  Exploring  both  simul¬ 
taneously  is  not  recommended  since  the  amount  of  data  would  prohibitive. 

The  nozzle  pressure  waveform  should  be  varied  both  in  form  and  maximum  am¬ 
plitude.  It  is  not  clear  what  influence  the  rise  in  pressure  has  upon  the  spray  charac¬ 
teristics.  This  study  could  only  vary  the  rise  while  varying  the  fall,  the  symmetric  case. 
Furthermore,  different  declining  pressure  profiles  are  now  required  to  quantify  and 
define  their  influence.  Changing  the  maximum  pressure  amplitude  may  change  the 
whole  issue.  The  maximum  pressure  used  in  this  study  is  between  low  pressure  gaso¬ 
line  injection  and  high  pressure  diesel  injection.  Therefore,  these  limits  must  be  ex¬ 
plored  for  practical  application. 

Resolving  a  particular  spray  in  higher  detail  is  important  to  the  validation  of  fu¬ 
ture  theories  and  models  regarding  intermittent  injections.  The  time-resolution  of  the 
Aerometric’s  PDPA  system  is  more  than  adequate  for  resolving  all  of  the  drops  pass¬ 
ing  the  sample  location  at  a  given  time.  However,  the  time-resolution  of  pressure 
should  be  improved  along,  and  instantaneous  injection  rate  should  be  monitored.  In 
addition,  more  sampling  locations  are  needed  to  completely  map  the  spray  and  its 
propagation,  verify  its  symmetry,  and  give  better  global  estimates  of  the  distributions 
for  the  injection  both  time-integrated  and  time-resolved. 
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Appendix  A 


Computer  Programs 

A.l  Analysis  of  Data  for  Local  Pressure  Peaks 

This  program  reads  through  the  pressure  data  (in  volts)  and  finds  local  pressure 
peaks.  Each  local  pressure  peak  represents  one  injection. 

Since  the  time  resolution  of  the  data  is  1  ms,  it  is  probable  that  the  actual  pres¬ 
sure  peak  occurs  between  samples.  An  interpolated  value  is  computed  from  the  re¬ 
corded  peak  and  its  two  neighboring  samples.  This  interpolated  peak  is  used  to  esti¬ 
mate  the  actual  peak.  A  second  order  polynomial  is  used. 

The  time  of  the  recorded  peak,  the  estimated  peak,  and  their  time  offset  from  the 
recorded  peak  are  then  saved  for  each  peak. 

The  conversion  from  volts  to  psi  is  done  by  the  graphing  software  (Harvard 
Graphics),  not  in  this  program.  The  conversion  is  a  simple,  linear  multiplier  (1  mv  = 
986  psi). 

PROGRAM  T1 

REAL  RED1 (2,20001), PPEAK , PEAKIN (3,500) 

LOGICAL  GOOD, GO, FOUND, STOP 

OPEN  (UNIT«1,FILE=* C:\DATA4\RAWP4.CON' ) 

OPEN  (UNIT-3,FILE='C: \DATA4\RED4\PEAK.DAT' ) 
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STOP* . TRUE . 

1*1 

BSUM-0 . 0 
ISUM*0 

WHILE  (STOP)  DO 

READ  (  1 ,  *  , IOSTAT-IOS ) RED1 (1,1), RED1 (2,1) 

XF  (XOS  .EQ.  0)  THEN 

IF  (RED1 ( 2 i X )  . 6E.  2)  THEN 

PRINT  *, 'LOOK' ,I,RED1(2,I) 

STOP 

END  XF 

IF  (RED1 (2,1) .LE.  .  t  THEN 
BSUM=BSUM+REDi v  i  ,  X  ) 

ISUM*ISUM+1 

END  IF 

ELSE 

STOP*. FALSE. 

END  IF 
I-I+l 
END  WHILE 
JMAX*I-2 

THRESH*. 1 

INTERPOLATE  PEAKS 
NUMPTS*0 
NUMPKS=1 
VPEAK=RED 1(2,1) 

IPEAK*RED1 (1,1) 

IND=1 

GOOD*. FALSE. 

DO  5  J*2 , JMAX 

XF  (RED1 ( 2 , J)  .GE.  THRESH)  THEN 
GOOD* . TRUE . 

IF  (RED1(2,J)  .GT.  VPEAK )  THEN 
VPEAK*R£D1 ( 2 , J ) 

IPEAK*RED1 ( 1 , J ) 

IND*J 

END  IF 

NUMPTS=NUMPTS+1 

ELSE 

IF  (VPEAK  .EQ.  R£D1(2, J-l) )  GOOD*. FALSE. 

IF  (VPEAK  .EQ.  RED1 ( 2 , J-NUMPTS ) )  GOOD*. FALSE. 
IF  ((GOOD)  .AND.  (NUMPTS  .GE.  5))  THEN 
F0-RED1 ( 2 , IND-1 ) 

C1*RED1 ( 2 , XND)-F0 
F1T1-F0+2*C1 

C2* ( RED1 ( 2 , IND+1 ) -F1T1 ) * . 5 
TPEAK-- ( 1+C1 /C2 ) * . 5 

PPEAK*F0+C1* (TPEAK+1 )+C2*TPEAK* (TPEAK+1 ) 
RECORDED  PEAK  INDEX  IN  RED1.DAT 

PEAKIN ( 1 , NUMPKS )-REDl ( 1 , IND ) 

OFFSET  BETWEEN  RECORDED  PEAK  TIME  AND  INTERPOLATED 
(-)  MEANS  INTERPOLATED  PEAK  COMES  FIRST 
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PEAKIN ( 2 , NUMPKS ) -TPEAK 


*  INTERPOLATED  PEAK 

PEAKIN ( 3 , NUMPKS ) -PPEAK 
NUMPKS-NUMPKS+1 
GOOD*. FALSE. 

END  IF 

VPEAK*REDI (  2 ,  J) 

TPEAK-RED1 ( 1 , J) 

IND-J 
NUMPTS*0 
END  IF 

5  CONTINUE 

WRITE  ( 3 , * ) NUMPKS -1 

DO  6  1*1, NUMPKS- 1 

WRITE  ( 3 , * ) PEAKIN (1,1), PEAKIN (2,1), PEAKIN (3,1) 

6  CONTINUE 

STOP 

END 

A.2  Tabulating  Drop  Size  and  Velocity  Data 

This  program  reads  through  the  list  of  local  pressure  peaks  created  by  program 
T1  (see  Appendix  A.1)  and  tabulates  drop  size  and  velocity  information  in  a  predefined 
time  window  (injection  event)  around  the  pressure  peak. 

The  correction  factors  for  Gausian  Errors  are  applied  to  the  final  distributions  in 
this  program.  A  list  of  correction  factors  is  computed  using  program  T7  (see  Appen¬ 
dix  A.4). 


PROGRAM  T2 

REAL  CORFAC (76) 

INTEGER  ISTORE( 120, 101 ) , ISTOR3 ( 120,101 ) , 

+  I DROPS ( 120 ) , ITOT( 76 ) , IDROPV ( 120 ) 

LOGICAL  GOOD, STOP 1, STOP 2 

OPEN  (UNIT*1,FILE*’C: \DATA4\RED4\PEAK.DAT' ) 

OPEN  (UNIT*2, FILE* 'C:\DATA4 \RED4\DISTR1.DAT' ) 

OPEN  (UNIT-3, FILE*’C: \PDPA\PDPDATA\6292\1\RUN04\PDPARAW.TXT'  ) 
OPEN  (UNIT-10, FILE-'C: \DATA4\RED4\VELDISTR.DAT' ) 

OPEN  (UNIT-4 , FILE- ' C : \DATA4 \RED4 \CORFACT.DAT ' ) 

DO  4  1-1,120 

DO  3  J-1,101 

IF  (I  .LE.  120)  ISTORE(I, J)-0 
IF  (I  .LE.  120)  ISTOR3 ( I , J)-0 
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3  CONTINUE 

IF  (I  .LE.  75)  ITOT(I)-0 
IDROPV(I)-0 

4  CONTINUE 

DO  5  1*1,75 

READ  ( 4 , * ) IDUM1 , CORFAC ( I ) 

5  CONTINUE 

READ  (3,*) 

READ  ( 1 , * ) NUMPKS , DUMMY 
NUMPKS*NUMPKS- 1 
STOP 1*. TRUE. 

I DROP* 1 
ICOUNT*l 

WHILE  ( STOP 1 )  DO 

READ  (1,*)RTIME, DUMMY, VOLT 
ITIME-RTIME 

IF  (ICOUNT  .BQ.  NUMPKS)  STOP 1*. FALSE. 

IF  ((VOLT  .GE.  .95)  .AND.  (VOLT  .LT.  1.05))  THEN 
ISTART-ITIME-10 
ISTOP*ITIME+50 
STOP 2*. TRUE. 

WHILE  ( STOP2 )  DO 

2  READ  ( 3 , * , IOSTAT*IOS )R1 ,R2 ,R3 ,R4 ,R5 

IT-INT(R5*1000 ) 

IF  ((IT  .LT.  ISTART)  .AND.  (IOS  .EQ.  0))  GOTO  2 
IF  ((IT  .GE.  ISTART)  .AND.  (IT  .LE.  ISTOP))  THEN 
ISIZE*INT(R1) 

IVEL*INT(R2 ) 

ITR-IT-ISTART+1 

ISTORE ( ITR, ISIZE )*ISTORE ( ITR, ISIZE )+l 
IF  ((IVEL  .LT.  76)  .AND.  (IVEL  .GT.  0)) 

+  ISTOR3 ( ITR, IVEL) -ISTOR3 ( ITR, IVEL )+l 

ISTORE (100, ITR )*ISTORE( 100, ITR )+l 
ISTORE ( 100,62 )*ISTORE( 100,62 )+l 

END  IF 

IF  ((IT  .GT.  ISTOP)  .OR.  (IOS  .NE.  0))  STOP2-. FALSE. 

IF  (IOS  .NE.  0)  STOP1-. FALSE. 

END  WHILE 
IDROP*IDROP+l 

END  IF 

ICOUNT*ICOUNT+l 
END  WHILE 

DO  15  1-1,61 
DO  14  J-1,75 

ISTORE ( I , J ) -ISTORE ( I , J ) *CORFAC ( J ) 

14  CONTINUE 

15  CONTINUE 

SUMI-C 
DO  7  1*1,61 

IDROPS ( I )*0 
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DO  6  J-1,75 

ITOT( J)-ITOT(J)+ISTOR£(I, J) 
DNORM-ISTORE ( I , J) / ( IDROP-1 . 0 ) 

*  WRITE  (2,*)J,DNORM 

IDROPS ( I ) -IDROPS ( I )+I STORE (I, J) 

6  CONTINUE 

*  WRITE  (2,*) 

6UMI-SUMI+IDR0PS ( I ) 

7  CONTINUE 

DO  11  1-1,61 

DO  10  J-1,75 

VNORM-ISTOR3(Z, J)/( IDROP-1. 0) 
IDROPV ( J ) « IDROPV ( J ) +ISTOR3 (I, J) 

*  WRITE  ( 10, •  )  J , VNORM 

10  CONTINUE 

11  CONTINUE 

ARNUM-0 
ARDEN* 0 
SANUM-0 
SADEN*0 
VENUM-0 
VEDEN=0 

DO  13  1*1,61 

DO  12  J=1 ,75 

ARNUM-ARNUM+ISTORE (2 ,J)*J 
ARDEN*ARDEN+ISTORE (I, J) 
SANUM—SANUM+ISTORE (I, J) * J**3 
SADEN*SADEN+ISTORE ( I , J ) * J**2 
VENUM=VENUM+ I STOR3 ( I , J ) * J 
VEDEN-VEDEN+ISTOR3 (I , J) 

12  CONTINUE 
ARMEAN-0 
SAMEAN-0 
VEMEAN-0 

IF  (ARDEN  .NE.  0.0)  ARMEAN-ARNUM /ARDEN 

IF  (SADEN  .NE.  0.0)  SAMEAN-SANUM/SADEN 

IF  (VEDEN  .NE.  0.0)  VEMEAN=VENUM/VEDEN 

ARNUM-0 

ARDEN -0 

SANUM*0 

S ADEN-0 

VENUM-0 

VEDEN “0 

DNORM-IDROPS ( I ) / ( IDROP-1 . 0 ) 

WRITE  ( 2 , * ) ARMEAN , SAMEAN , DNORM 
WRITE  ( 10 , * ) I , VEMEAN 
13  CONTINUE 

WRITE  ( 2 , * ) SUMI 

DO  20  J-1,75 

DNORM-ITOT(J)/( IDROP-1. 0) 
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VN0RM-IDR0PV(J)/(IDR0P-1. ) 

WRITE  (2 ,  * ) J,DNORM 
*  WRITE  ( 10 , * )  J , VNORM 

20  CONTINUE 

STOP 

END 

A3  Analysis  of  Pressure  Waveforms 

This  program  calulates  average  pressure  waveforms  for  each  of  the  pressure  wave¬ 
forms  used.  In  addition,  a  standard  deviation  is  computed  for  each  time  point  in  a  give 
pressure  waveform.  The  standard  deviation  measures  how  well  any  individual  prussure 
curve  might  compare  to  the  average  at  a  give  time  point. 


PROGRAM  T6 

REAL  STORE ( 120 ) , STORE1 ( 120 ) 

LOGICAL  GOOD , STOP1 , STOP2 

OPEN  (UNIT-1, FILE-'C: \DATA4\RED10\PEAK.DAT' ) 

OPEN  (UNIT-3, FILE- 'C:\DATA4\RAWP10. CON ' ) 

OPEN  (UNIT-2, FILE- 'C:\DATA4\RED10\PRESS.DAT' ) 

OPEN  (UNIT-4, FILE«'C: \DATA4\RED10\PRSTAT.DAT' ) 

READ  (2,*) 

DO  4  1-1,41 

READ  ( 2 , * ) DUM , STORE ( I ) 

STORE 1 ( I )-0 . 0 
4  CONTINUE 

READ  (3,*) 

READ  ( 1 , * )NUMPKS 

READ  ( 1, *) RTIME, AD J, VOLT 

NUMPKS-NUMPRS-1 

STOP1-.TRUE. 

IDROP-1 

ICOUNT-1 

WHILE  (STOP1)  DO 

READ  ( 1 , * ) RTIME , AD J , VOLT 
ITIME-RTIME 

IF  (ICOUNT  .EQ.  NUMJKS)  STOP 1-. FALSE . 

IF  ((VOLT  .GE.  .95)  .AND.  (VOLT  .LT.  1.05))  THEN 
WRITE  (4,*) RTIME, VOLT 
ISTART-ITIME-10 
ISTOP-ITIME+30 
8TOP2-.TRUE. 

READ  (3,*) 

IF  (ADJ  .LE.  0.)  THEN 
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WOLTP-O . 

WHILE  ( ST0P2 )  DO 

2  READ  ( 3 ,  * , IOSTAT*IOS  )TIME,  WOLT 

IT=TIME 

IF  ((IT  .LT.  I  START)  .AND.  (IOS  .BQ.  0))  GOTO  2 
IF  ((IT  .<Z.  I  START)  .AND.  (IT  .IE.  ISTOP) )  THEN 
I TR* I T- I START+ 1 

AD  VOLT*  ( WOLT-WOLTP )  *  ( 1 .  -ADJ )  +WOLTP 
IF  (IT  .BQ.  ISTART+10)  ADVOLT-VOLT 
siazi(m)^iazi(m)^AiwxT-SKrE(iiR)  >**2 

END  IF 

IF  ((IT  .GT.  ISTOP)  .CR.  (IOS  .NE.  0))  STCP2».PALSE. 
IF  (IOS  .NE.  0)  STOP1*. FALSE. 

WOLTP-WOLT 
END  WHILE 

ELSE 

READ  (3,*) TIME, WOLT 
WHILE  (STOP2 )  DO 

6  READ  ( 3 , * , IOSTAT*IOS  )TIMEA,  WOLTA 

IT*TIME 

IF  ((IT  .LT.  ISTART)  .AND.  (IOS  .EQ.  0))  THEN 
TIME=TIMEA 
GOTO  6 

END  IF 

IF  ((IT  .GE.  ISTART)  .AND.  (IT  .LE.  ISTOP))  THEN 
ITR-IT-ISTART+1 

ADVOLT*  (  WOLTA- WOLT )  *  ( ADJ )  +WOLT 
IF  (IT  .EQ.  ISTART+10)  ADVOLT=VOLT 
STORE  1  ( ITR ) “STORE  1  ( ITR )+< ADVOLT- STORE (ITR) )**2 

END  IF 

IF  ((IT  .GT.  ISTOP)  .OR.  (IOS  .NE.  0))  STOP2=.FALSE. 

IF  (IOS  .NE.  0)  STOP1*. FALSE. 

WOLT=WOLTA 
TIME=TIMEA 
END  WHILE 

END  IF 

IDR0P=IDR0P+1 

END  IF 

ICOUNT«ICOUNT+l 
END  WHILE 

WRITE  ( 4 , * )IDROP-l 
DO  7  1*1,41 

DNORM-SQRT( STORE1 ( I ) / ( IDROP-1 . 0 ) ) 

WRITE  ( 4 , * ) I , DNORM 

7  CONTINUE 

STOP 

END 

A.4  Calculation  of  Correction  Factors 

This  program  reads  the  corrected  and  uncorrected  distributions  created  by  the 
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PDPA  software  and  calculates  the  correction  factors  used.  The  correction  factors  are 
then  used  to  create  the  list  of  correction  factors  for  the  different  size  groups  in  the  cus¬ 
tom  distributions. 

PROGRAM  CORREC 

*  CALCULATION  OF  CORRECTION  FACTORS  TO  BE  USED  WITH  CUSTOM 

*  DISTRIBUTIONS. 

REAL  STORE ( 3,76) 

INTEGER  CORVAL 

OPEN  (UNIT*2 ,FILE= ' C: \PDPA\PDPDATA\6292\1\RUN23\PDPADATA.TXT'  ) 

OPEN  ( UNIT=1 , FILE* ’ C  s \ DATA 4 \RED2  3 \CORFACT . DAT ' ) 


DO  4  1*1,3 
DO  3  J-1,75 

STORE(I, J)*0 

3  CONTINUE 

4  CONTINUE 

DO  5  1=1,72 

READ  (2,*) 

5  CONTINUE 

DO  6  1=1,50 

READ  ( 2 , * ) IDUM1 , STORE (1,1), STORE (2,1) 

6  CONTINUE 

READ  (2,*) 

READ  (2,*) 

DO  7  1=1,50 

READ  ( 2, *)IDUM1,DUM2, CORVAL 

IF  ( STORE (2,1)  .NE.  0.)  STORE ( 2 , 1 )*CORVAL/ STORE ( 2 , 1 ) 

7  CONTINUE 

J=1 

STORE ( 3 , 1 )*STORE ( 2 , J) 

WRITE  ( 1 , * ) J , STORE ( 3 , J ) 

DO  8  1=2,75 

AVBETW* ( STORE ( 1 , J)+STORE(l, J+l) )/2. 

IF  (I  .LE.  AVBETW)  J*J 
IF  (I  .GT.  AVBETW)  J«J+1 
STORE (3,1) “STORE ( 2 , J ) 

WRITE  ( 1 , * ) I , STORE (3,1) 

8  CONTINUE 


STOP 

END 
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26797-EG 


S.P.  Lin  Clarkson  University 


C  APPfcNpoc  f  ) 


"Mechanism  of  Intermittent  Atomization" 


This  research  project  seeks  a  complete  knowledge  of  the  physical  mechanism  of  the 
formation  of  intermittent  fuel  sprays.  The  objective  is  to  exploit  this  knowledge  to  enhance  the 
efficiency  of  internal  combustion  engines,  and  to  reduce  the  air  pollution  caused  by  inefficient 
combustion.  The  theoretical  results  (S.P.  Lin  10th  ARO  Engine/Fuel  Workshop,  Proceedings  pp. 
120-131,  Madison,  Wisconsin,  June  2-3,  1992)  show  that  the  droplet  formation  from  the  leading 
edge  of  the  liquid  fuel  jet  is  caused  by  the  Taylor-Lamb  instability  during  the  decerelation  phase 
of  the  intermittent  spray.  The  mechanism  of  droplet  formation  from  the  rest  of  the  liquid-gas 
interface  remains  the  same  as  that  for  the  case  of  stationary  atomization  (S.P.  Lin  and  Z.W.  Lian, 
pp.  771-773,  Physics  of  Fluids  A,  1993)  Since  the  leading  edge  interfacial  area  occupies  only 
a  small  part  of  the  total  interfacial  area,  the  atomization  of  intermittent  sprays  is  still  caused  by 
the  pressure  fluctuation  at  the  interface.  However,  this  does  not  imply  that  the  size  distributions 
for  stationary  and  intermittent  sprays  are  the  same.  In  fact  the  decerelation  of  the  liquid  jet  tends 
to  reduce  the  droplet  diameter  at  the  same  pressure.  To  compare  with  the  theory,  the  size  and 
velocity  distributions  at  various  locations  in  a  spray,  as  indicated  in  Fig.  2.41,  are  measured  with 
an  Aerometric  Phase  Dopier  Particle  Analyzer.  The  figure  numbers  refer  to  the  figure  numbers 
of  the  M.S.  thesis  of  V.  Cook  (V.F.  Cook,  M.S.  Thesis,  Clarkson  University,  1993).  Fig.  3.3.3 
gives  the  Sauter  mean  diameters  at  the  location  B-A  as  function  of  time  for  3  different  sprays 
with  different  forms  of  pressure  pulse  at  the  nozzle  inlet.  All  3  sprays  has  the  same  peak 
pressure.  One  spray  has  a  symmetric  pressure  pulse  of  period  14  ms.  The  second  spray 
designated  with  steep  TE  has  a  sudden  pressure  reduction  after  the  pressure  peak.  The  third 
spray  designated  with  gentle  TE  has  a  gentle  pressure  reduction  after  the  pressure  peak.  It  can 
be  seen  that  the  Sauter  mean  diameter  is  the  smallest  on  the  average  for  the  case  of  steep  TE. 
This  has  profound  practical  implications.  The  corresponding  size  distributions  are  given  in  Fig. 
3.1.1.  Comparing  the  number  of  drops  from  the  gentle -TE-spray  and  the  steep-TE-spray  in  this 
figure,  it  is  seen  that  the  former  spray  produces  more  droplets  than  the  latter  one.  Moreover  the 
former  produces  relatively  more  smaller  droplets  than  the  latter.  This  indicates  that  the 
mechanisms  of  intermittent  and  stationary  spray  formation  remain  essentially  the  same.  This  is 
because  the  gentle  TE  spray  has  a  larger  velocity  over  the  entire  period  of  deceleration,  and  it 
is  known  (Z.W.  Zhou  and  S.P.  Lin,  J.  of  Propulsion  and  Power,  8,  pp.  736-740,  1992)  that 
smaller  droplets  can  be  produced  with  larger  jet  velocities  in  stationary  jets.  To  support  this 
argument,  the  size  distributions  of  an  intermittent  sprays  at  different  times  during  the  period  of 
pulsation  were  measured  and  compared  with  the  stationary  spray  data  for  different  steady  nozzle 
pressures.  The  latter  data  are  given  in  Fig.  3.4.1a.  The  data  for  the  intermittent  spray  are  not 
given  here,  because  of  the  lack  of  space.  The  comparison  reveals  that  the  shape  of  size 
distributions  curve  of  an  intermittent  spray  changes  according  to  the  change  in  pressure  in  the 
same  manner  as  that  for  the  stationary  sprays.  It  can  be  shown  theoretically  that  intermittent 
sprays,  as  far  as  the  mechanism  of  droplet  formation  is  concerned,  remain  quasi-steady  unless 
the  period  of  pulsation  is  less  than  100-th  of  an  ms  for  diesel  fuel.  However,  finer  sprays  can 
be  produced  w'ith  a  more  sudden  reduction  in  the  nozzle  inlet  pressure  after  its  peak  during  the 
injection  period. 


Nozzle 


X  (radial  direction) 


•  B-B  (1 32cm,  8.4cro) 

•I 

B-A  (132cra,23.9cra) 


A-A  (032cra,23.9cm) 


•••  *D-A  (232cm,  23.9cm) 

C-A  (0.82cra,23.9cra) 

Z  (axial  direction) 


Figure  2.4.1  Spray  sampling  locations  downstream  of  the  nozzle . 
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Figure  3.1.1  Time-integrated  drop  size  distributions 
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Figure  33.3  Time  evolution  of  the  Sauter  mean  diimeter  1 
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Figure  3.4.1  Continuous  789,  888,  and  986  psi  nozzle  pressure  spray  drop  size  (a)  and  drop  velocity 
(b)  distributions 


The  absolute  and  convective  instability  of  a  viscous  liquid  jet  emanating  into  a  viscous 
gas  in  a  vertical  pipe  is  analyzed  in  a  parameter  space  spanned  by  the  Reynolds  number,  the 
Froude  number,  the  Weber  number,  the  viscosity  ratio,  the  density  ratio,  and  the  diameter  ratio. 
The  numerical  results  of  the  analysis  are  used  to  demonstrate  that  reduction  in  gravity  tends  to 
enhance  the  Rayleigh  mode  of  convective  instability  which  leads  to  the  breakup  of  a  liquid  jet 
into  drops  of  diameters  comparable  with  the  jet  diameter.  On  the  contrary,  the  Taylor  mode  of 
convective  instability  which  leads  to  atomization  is  retarded  at  reduced  gravity.  The  Rayleigh 
mode  becomes  absolutely  unstable  when  the  Reynolds  number  exceeds  a  critical  value  for  a 
given  set  of  the  rest  of  the  relevant  parameters.  The  domain  of  absolute  instability  is 
significantly  enlarged  when  the  effect  of  gas  viscosity  is  not  neglected. 

A  new  theory  is  developed  for  an  intermittently  atomizing  jet  emanating  from  a  nozzle 
into  a  combustion  chamber  of  finite  dimension.  The  velocity  and  pressure  fields  in  the  nozzle, 
jet,  and  the  chamber  are  described  by  Eulerian  description.  The  description  is  translated  into 
Lagrangian  particle  description  at  the  liquid-gas  interface.  Thus,  the  computation  domain  is 
gridless  but  is  rapidly  increasing  and  decreasing  w'ith  time.  So  far  as  the  writer  is  aware,  this  is 
the  first  theory  which  is  capable  of  revealing  the  atomization  process  in  response  to  the  pressure 
input.  There  are  two  major  phases  in  intermittent  atomization:  acceleration  and  deceleration 
phases.  During  the  acceleration  phase,  atomization  initiates  from  the  size  of  a  jet.  During  the 
deceleration  phase  more  droplets  are  formed  from  the  penetrating  front  of  the  jet.  The  droplets 
formed  during  the  former  phase  is  smaller  than  those  formed  in  the  latter  phase.  This  is 
demonstrated  in  the  attached  figure.  The  time  evolution  of  interfacial  positions  of  2  dense  sprays 
for  2  different  pressure  pulses  introduced  at  the  nozzle  inlet  are  delineated  in  this  figure. 

I  visited  Dr.  W.  Bryzik  of  TACOM  on  21  September  1992.  The  visit  made  me  aware  of 
a  fundamental  problem  in  atomization  which  is  even  more  basic  than  the  intermittency  of  sprays. 


